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SUMMARY 


This report summarizes the results of a 12-month study of real-time 
computer simulation requirements for the Rotor System Research Aircraft 
(RSRA). The report develops scaled equations suitable for simulation of the 
rotor system. These equations form the basis for consideration of both digital 
and hybrid mechanization for real time simulation. For all-digital simulation 
estimates of the required speed in terms of equivalent operations per second are 
developed based on the complexity of the equations and the required integration 
frame rates. For both conventional hybrid simulation and hybrid simulation 
using time -shared analog elements the amount of required equipment is 
estimated along with a consideration of the dynamic errors. It is concluded 
that conventional hybrid mechanization using analog simulation of those rotor 
equations which involve rotor -spin frequencies (this constitutes the bulk of the 
equations) requires too much analog equipment. Hybrid simulation using time- 
sharing techniques for the analog elements appears possible with a reasonable 
amount of analog equipment. All-digital simulation with affordable general- 
purpose computers is not possible because of speed limitations, but specially- 
configured digital computers do have the required speed and constitute the 
recommended approach. 
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1. INTRODUCTION 


This report presents the results of a 12 -month study of real-time 
computer simulation requirements for the Rotor System Research Aircraft 
(RSRA). The original math model to be used as a basis for the study included 
only the equations for the rotor system and was furnished by the Simulation 
Sciences Division of the NASA Ames Research Center^ These equations 
were originally developed by Sikorsky Aircraft and have also been used 
extensively by the NASA Langley Research Center. At both Langley and Ames 
major difficulties have arisen in trying to simulate the RSRA in real time, 
largely because of the complexity of the rotor model. This motivated the 
present study, which originally had the objective of developing a hybrid 
computer mechanization of the rotor equations in order to permit adequate 
real-time simulation. 

In order to better understand the RSRA equations, and in particular 
the rotor system equations, the simplifications in the equations and the 
frequency content in the various state variables, the equations were rederived 
from first principles (see Appendix A) and rewritten in scaled, dimensionless 
form (Appendix B). In addition to the rotor -system equations the six-degree- 
of -freedom rigid aircraft equations have also been included with the exception 
of the aerodynamic terms. Thus the only aerodynamic forces considered in 
the model are the aerodynamic forces acting on the rotor blade segments. 
These forces, however, are predominant and also the complexity associated 

^Mackie, D. Brian and Thomas S. Alderete, "A Real-Time Dual Processor 
Simulation of the Rotor System Research Aircraft, ” NASA TN-D-8328, 
January, 1977. 
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with calculating them and the resulting blade motion and forces is predominant 
in the computer simulation of the RSRA. For this reason the RSRA model 
considered in this report is felt to be reasonable as a basis for comparing 
different computer mechanizations. The additional computation required for 
the fuselage, tail, and wing aerodynamic forces should not change substantially 
the conclusions of the study. This is especially true when one considers that 
the rotor blade equations in their entirety are not only more complex but also 
contain much higher frequencies than the aircraft equations. 

In the rotor model each blade is represented as a rigid body with two 
degrees of freedom (lagging and flapping). The model can be altered to allow 
flexible blades by using normal coordinates for the elastic degrees of freedom, 
but this would further increase the equation complexity and high-frequency 
content of the problem. It is probably not necessary for simulation associated 
with either man-in-the-loop or control- system studies. 

In Chapter 2 the scaled equations in Appendix B are used as a basis 
for a total count of various computational operations. This count, in turn, 
allows requirement estimates for various types of computer mechanization 
of the equations. Chapter 3 discusses the requirements for all-digital 
mechanization, Chapter 4 for conventional hybrid mechanization and Chapter 5 
for hybrid mechanization using a time-shared analog system for computing 
blade forces. Finally, Chapter 6 provides summary recommendations 
as a result of the study. 
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2. COMPUTATIONAL OPERATIONS REQUIRED TO IMPLEMENT 

THE RSRA MODEL 

2. 1 Basic Model 

The math model for the RSRA is developed in Appendix A. Figure 2. 1 
in Section A. 2 in the appendix shows a block diagram of the equations in vector 
form. The six-degree -of-freedom rigid aircraft equations are shown in block 1 
of the figure. The equations for the two degrees of freedom (lagging and flapping) 
for each of the rotor blades are shown in blocks 9 and 12. Blocks 2 through 8 
represent the equations for transforming the angular and translational velocity 
and acceleration components of the aircraft to velocity and acceleration 
components along each rotor blade. Block 10 represents the computation of 
aerodynamic forces (lift and drag) on each blade segment and block 11 represents 
the calculation of the resulting moments on each blade, as needed for solving 
the blade equations of motion. Finally, blocks 13 through 15 represent the 
conversion of individual blade forces to overall force and moment components 
acting on the rigid aircraft. 

In Appendix A the scalar equations in each of the 15 blocks in Figure 2. 1 
are developed in the 15 respective subsections of A. 2. In Appendix B these 
equations are rewritten in dimensionless form as needed for analog /hybrid or 
fixed-point digital solution. It is the equations of Appendix B which are used as 
the basis for determining real-time computer requirements. Certain 
modifications, such as subtracting off the principal centrifugal acceleration 
term in computing blade accelerations, are made in the equations in Appendix B 
to improve scaling and avoid taking small differences of large quantities. 

These modifications appear necessary in order to use the RSRA rotor model 
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in analog or 16-bit fixed-point digital simulation. They are probably not 
needed for floating-point or 24 bit fixed-point digital computation. The 
complexity added to the equations with these modifications (the equations are 
still exact) is minimal. 

The equations in Appendix B do not have any terms dropped for simpli- 
fication purposes. This could be done in a number of cases (e. g. , dropping 
terms involving the product of sines of two small angles) but the resulting 
simplifications in the equations would not be very substantial. In general the 
equations are similar to the ARC equations^ except in the approach for computing 
and /3^, the time-rate-of-change of the lagging and flapping angles, 

respectively, for the blade. In the ARC equations 6. and j3. are computed 

J 

by approximate equations of motion, from which 6^., 6^, and are computed 
assuming a periodic fourier series form for the solution. In Appendices A 
and B the exact blade equations are integrated to compute blade angular 
velocity components, which are then subtracted from rotating shaft-axis 
angular velocity components to compute 6^. and )3. respectively. These in turn 
are integrated to obtain 6^ and The equation complexity is comparable with 
the ARC equations. Actual computer simulation will be needed to determine 
which method gives better results. 

Table 2, 1 summarizes the operations count involved in each of the 
computational blocks. The operations shown are needed to compute the right- 
hand side of all equations, including those with first time derivitives on the 

^Mackie and Alderete, NASA-TN- D-8328, op. cit. 
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left side. Thus the additional operations needed to implement the numerical 
integration formulas are not included in the count. This is because those 
operations are simple and qiiite small in number compared with the operations 
required to complete each pass through the calculation of all the state -variable 
derivitives, which is what Table 2. 1 does represent. The table shows the 
count based on N rotor blades with s segments per blade. Where the equations 
must be solved for each blade, the operations count is multiplied by N; where 
the equations must be solved for each segment of each blade, the operations 
count is multiplied by Ns. 

It should be noted that there are some alternate ways of mechanizing some 

of the equation terms. For example, in Eq. (12. ]) the term (1 - cos is 

counted as a one-variable function, whereas it might be acceptable to compute 

2 

the term approximately as 0.5^^ . Similarly, the term (1 - cos 6^ cos 

in Eq. (13. 2) is counted as a two-variable function, whereas it might be accurate 

2 2 

enough to compute it approximately as (6^ + /3^ )/2. Also, in Eq. (10. 14), 
Appendix A, the equation for segment angle of attack, covinted as 

a two -variable function, whereas it can alternatively be mechanized as a 
combination of boolean and algebraic computations^. 

In the count of multiple /divides in Table 2.1a separate column is 
shown for digital and analog mechanization. This is because multiplications of 
variables by constants require multiply operations in digital computation, 
whereas they can be implemented using coefficient devices in analog computation. 


Mackie and Alderete, NASA-TN D-8328, op. cit. 
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Table 2. 1 Summary of Operations Count for Solving the RSRA Equations for 
N Blades, s Segments per Blade 


Block #, 

Eq. #in #Mult/Div's 

# # Trig 

#1 Var 

#2 Var 

# Sq. 

Fig 2. 1, Ap.A 

Ap. B Digital Analog 

Adds Fncts 

Fnct's 

Fnct's 

Roots 

1. 

(1.7-1, 9) 

14 

8 

9 

4 





(1. 12-1. 14) 

9 

3 

3 






(1. 16-1. 18) 

6 

6 

3 






(1, 19-1. 21) 

21 

21 

10 

2 




2. 

(2,2-2. 4) 

8 

0 

5 






(2.5-2. 7) 

4N 

4N 

2N 

2N 





(2. 8-2. 9) 

7N 

7N 

3N 

4N 




3. 

(3. 1-3. 3) 

8 

0 

5 






(3.4-3. 5) 

2N 

2N 

N 





4. 

(4. 3-4. 5) 

6 

0 

3 






(4. 6-4. 8) 

11 

0 

8 





5. 

(5. 3-5. 5) 

15 

6 

12 






(5. 7-5. 9) 

8 

0 

5 





6. 

(6. 1-6, 3) 

7N 

4N 

5N 





7. 

(7. 1-7. 3) 

4N 

4N 

2N 






(7.4-7, 6) 

12N 

6N 

9N 





8. 

(8. 1-8. 3) 

ION 

ION 

5N 






(8. 6-8. 8) 

6 Ns 

0 

3Ns 





9. 

(9. 1-9.3) 

12N 

ION 

7N 






(9. 5-9. 6) 

9N 

3N 

6N 





10. 

(10. 7-10. 9) 

llNs 

lONs 

3 Ns 






(10.3,4,ApA) 

4Ns 

4Ns 

Ns 




Ns 


(10.5) 

Ns 

Ns 

Ns 




Ns 


(10. 23,ApA) 

6N 

N 

N(5+s) 

2N 





(10. 12,ApA.) 

6Ns 

6Ns 

2 Ns 

2Ns 





(10. 15,ApA.) 

2 Ns 

Ns 

Ns 




Ns 


(10. 14, 16,17, ApA.) 





3Ns 


11, 

(11, 1-11. 2) N(3)2s) 0 

N(242s) 

N 



12. 

(12. 1-12.2) 

4N 

2N 

3N 


N 



13. 

(13. 1-13.3) 18N 

8N 

7N 



N 



(13. 6-13. 8) 

3N 

0 

3 Ns 





14. 

(14. 1-14.3) ION 

ION 

7N 






(14.4-14.6) 

6N 

6N 

4N 






(14. 8-14.10) 

N 

0 

3(N-1) 






(1411-14.13) 

8 

0 

5 





15. 

(15. 2-15.3) 

3N 

N 

2N 






(15. 5-15. 7) 

2N 

2N 

0 






(15, 8-15.10) 

3N 

0 

3(N-1) 






(15.11-15.13) 

8 

0 

5 






(15.15-15.17) 

9 

0 

6 





Total Operations Count 

131 

44 

73 

6 






+ 

+ 

+ 

+ 






126N 

SON 

75 N 

8N 

2N 

N 




+ 

+ 

+ 

+ 


+ 




32Ns 22Ns 

18Ns 

2Ns 


3Ns 

3Ns 
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At the end of Table 2. 1 is shown a total count for each of the operation 


categories. Where a count in multiplied by N, the number of operations 
depends on the number of rotor blades, N. It also means that those equations 
involve the rotating blade frequency typically the order of 30 radians per 
second, and N can range up to 6. Where a count is multiplied by Ns the nixmber 
of operations depends not only on the blade count N, but also the number of 
segments s taken per blade. A typical value for s is 5, so that Ns might be 
6x5 = 30. Hence those operations multiplied by Ns can become very 
numerous indeed. They also clearly involve the rotor frequency 

The total operations count for each column in Table 2. 1 is used in the 
next three sections to estimate the requirements for real time digital and 
hybrid simulation. It shotald be noted that the equations in Appendices A and 
B and hence Table 2. 1 do not include the aerodynamic forces due to aircraft 
fuselage, wing and tail, nor do they include rotor down wash effects. It is 
believed that this omission does not affect substantially the conclusions in the 
next three sections. 
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3. ALL-DIGITAL SIMULATION 


3. 1 Choice of Integration Method 

The choice of numerical integration method for implementing an 
all-digital simulation of the RSRA is important in establishing the re- 
quired integration frame rates and the corresponding estimate of digital 
operations per second based on Table 2,1. For non real-time helicopter 

simulation, studies have shown that the 4th-order Runge-Kutta method is a 
2 

good choice . Although RK-4 gives good accuracy for fairly large integra- 
tion step size and also has good stability characteristics, it does require 
four passes through the equations (four derivative evaluations ) per integra- 
tion step. It also requires input data points one -half and one step ahead 
of the current point in time for evaluation of the required derivatives. 

This is a problem in real time simulation with external inputs \mless ex- 
trapolation algorithms are used on the inputs. 

If dynamic accuracies of the order of one percent are adequate for 
RSRA flight simulation, which we believe they are, then 2nd-order Rimge 
Kutta (RK-2) may be preferable over fourth order (RK-4). It still suffers 
from the real-time data input problem noted above, however. 

The data input problem is eliminated by using a mtiltistep predictor 
algorithm such as Adams -Bashforth. Such numerical integration algorithms 
have the added advantage of requiring only one derivative evaluation per 

2. Davis, J. M. , R.L. Bennett, and B. L. Blankenship, "Rotor craft 
Flight Simixlation with Aeroelastic Rotor and Improve Aerodynamic 
Representation", USAAMRDL-TR-74-10A, Fort Eustis, Va. , Jime, 1974. 
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integration step, the other derivatives needed for the algorithm having 
been already evaluated in previous steps. Also, for dynamic accuracy 
in the area of one percent, 2nd order Adams-Bashforth may be a good 
compromise and has indeed been widely used in digital flight simulation. 

In lieu of actually trying various integration schemes on the full 
RSRA equations, the most effect analysis technique is to consider lin- 
earized equations about some trim or steady -state condition. In this case 
the overall linearized set of equations will have n characteristic roots, 
i. e, the system transfer fxxnction will have n poles. With the use of a 
partial fraction expansion the system can then be decomposed into separate 
paralleled first and second order linear systems. The errors in transfer 
functions and/or characteristic roots due to any integration algorithm can 
be studied by considering the errors in the individual first and second 
order systems due to each algorithm. Then by superposition the overall 
dynamic error can be estimated. 

3 

For example, with z -transform analysis methods we can show that 
using RK-2 integration, the damping ratio and fractional frequency errors 
in simulating an under damped second -order system are the following: 

^ 7174 ^ 

i - e = -fr) (aT^ - c^ - e + I ) (aT)^ (3.1) 

A 

f - 1 = (^-|^^) (aT)^ , aT «1 (3.2) 

3. Gilbert, E, G. , "Dynamic Error Analysis of Digital and Combined 
Analog -Digital Computer Systems", SIMULATION, Vol. 6, No. 4, April, 
1966, pp* 241-257. 
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y 


Here ^ is the exact damping ratio and ^ the damping ratio resulting from 

A 

use of RK-2. Similarly, co is the exact root frequency and co the frequency 


when using RK-2. T is the integration interval and a the undamped natural 


frequency of the system. For given accuracy requirements in ^ and co 


Eqs. (3.1) and (3. 2) can be used to compute aT and hence the number of 


frames (integration steps) per cycle of the transient oscillation exhibited 


by the second-order system. 


As an illustration, let aT = 0.2, This corresponds to 5 integration 


steps per radian or 31.4 steps per cycle. Let ^ = 0,1 and assume we use 

A 

RK-2. Then from Eqs. (3. 1) and (3. 2) we find that ^ ^ = 0. 001 (damping 

A 

CO 

ratio is . 101 instead of . 1) and “ -1 = . 0064 (frequency is 0. 64% high). 


The counterparts of Eqs. (3.1) and (3.2) are more complex for 

3 

RK-4. For ^ = 0 the formulas are given by 


i=0 I - e = e = I (aT)^ (3.3) 

A 

^ -l = -^(aT)^ aT«l (3.4) 


For the case considered above for RK-2, namely ^ , 1, numerical solu- 


tion of the z-transform equations shows that with RK-4 using aT = 0.4 

A 

(2, 5 integration steps per radian or 15.7 per cycle), ^ - ^ = 0, 00059, 

A 

CO /co-1 = -0. 0036 (3. 6%). In this case RK-2 and RK-4 give comparable 


accuracy when the step-size for RK-4 is twice as big, which actually 


means the same number of passes (derivative evaluations) per integration 

4 2 

step. Since RK-4 errors depend on (aT) compared with (aT) for RK-2, 


3, Gilbert, op^ cit 
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higher accuracy simulation requiring smaller step sizes would favor use 


of RK-4. 

For End-order Adams -Bashf or th (AB-E) the following formula 
counterparts to Eqs. (3. l)and (3, E) can be derived: 


e - e 


43 

3E 


169 a 
96 ^ 


(aT)" 


A J- E09 .2 , 169 .4 

CO 2 1^ " 96 ^ 96 ^ 

o> ^ , ^2 


(aTr, aT « 1 




(3.5) 


(3.6) 


Using AB-E in simulating a second order system with damping ratio 4 = 0. 1 

and aT = 0.1 (62. 8 integration steps per cycle) we find from numerical 

A A 


solution of the z -transform equations that 4 -4=0- 00101 and cj /co -1 = 

0, 00412 (0.4%). The damping ratio error is the same as that for RK-2 
with aT = 0.2 (31.4 steps /cycle) and the frequency error is about 2/3 that 
of RK-2, Both schemes under these conditions require the same number 
of passes (derivative evaluations) per integration step. This is because 
RK-2 requires two passes per step versus one per step for AB-2. 

On the basis of the above type of analysis we conclude that for 
dynamic accuracy of the order of one percent End order Adams -Bashforth 
with aT somewhere between 0. 1 and 0. 2 will be satisfactory, i. e. , say 
~ 50 integration steps per cycle of the highest important problem frequency. 
One difficulty with Adams -Bashforth is the startup problem, but for an 
ongoing real-time simijlation this shoxaldn* t result in significant errors. 
3rd-order Adams -Bashforth should probably also be considered and 
analyzed. 
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3. 2. Required Integration Frame Rates 


If we settle on 2nd-order Adams Bashforth as a proposed integra- 
tion scheme with 50 steps per cycle of highest problem frequency, then 
the next consideration in real-time digital simulation is an estimate of the 
highest problem frequencies in the RSRA. For the rigid aircraft itself 
1 hertz is a reasonable estimate. When non -elastic rotor blades are 
considered, then the maximum rotor speed, roughly 30 radians per second 
or about 5 hertz, is a dominant high frequency. This suggests a minimum 
of 5 X 50 = 2 50 integration steps per second (4 millisecond step size) 
using AB-2 wherever the rotor -blades are involved in the equations. In- 
cidentally, 50 integration steps per cycle corresponds to one step every 

7.2 degrees of azimuthal advance for the rotor, a result not inconsistent 

4 

withHouck*s conclusions , His studies showed that somewhat larger 
steps can be taken, but he was using a higher-order integration algorithm 
(RK-4, we believe). There has been some speculation that the integration 
frame rate shoiild be based on a highest problem frequency equal to the 
rotor speed times the number of blades. Nevertheless, we will assume 
a 4 millisecond integration step using 2nd order Adams -Bashforth (one 
pass per step) to analyze the real-time digital computer requirements. 

3. 3 Required Digital Operations per Second for the Real-Time RSRA 

Simulation 

Let us assume 6 rotor blades with 5 segments per blade as 

4. Houck., J. A, , Computational Aspects of Real-Time Simulation of 
Rotary-Wing Aircraft”, Masters Thesis, George Washington University, 
May, 1976. 
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representative of a worst case RSRA simiilation. Then N = 6 cind s = 5 in 
Table Z.l, which results in the breakdown of digital operations per pass 
shown in Table 3.1. Clearly an overwhelming fraction of the total opera- 
tions are dependent on the blade number, N, and hence involve the rotor 
spin frequency of 30 rad /sec. Very little would be saved by solving the 
non-N dependent equations at a lower frame rate. 

From Table 3.1 rough estimates of actual frame times required 
for one pass through the equations and hence one integration step (assum- 
ing AB-2) can be made for specific digital computers. The actual number 
of individual digital operations (load, store, add, miiltiply, sign-test, 

etc. ) for each of the categories in Table 3.1 can be estimated, assuming 

5 

a conventional general purpose computer with no pipelining . For ex- 
ample , to compute a function of two variables takes 34 operations (10 adds, 

6 multiply /divide s , 13 loads, 1 store, and a nominal 4 sign tests). Since 
the multiply /divide time on many general-purpose computers takes roughly 
4 times the time required for the other operations, we use the concept of 
normalized operations by weighting all multiply /divides by a factor of 4. 
Using this concept the above two-variable function computation takes 46 
normalized operations per pass. Extending this concept to the total coiint 
in Table 3.1 yields 23,603 normalized digital operations required per pass 
through the RSRA equations. Since typical general purpose computers in use 

5. Howe, R. M, , and J. Paul Landauer , *'A Quantitative Method of Speed 
Comparison between Analog/Hybrid and Digital Computers", IEEE Com- 
puter Magazine, Vol. 9, No. 7, July, 1976, pp. 31-36. 
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Table 3.1 Digital Operations per Derivative Evaluation for N = 6 
s = 5 Segments per Blade 


# # 

M-H/Div^s Adds 

Operations 

Independent 131 73 

of N or s 


#Trig #1 Var #2 Var 
Fnct* s Fnct* s Fnct* s 

6 


Operations 

dependent 756 450 48 12 6 

on N = 6 


Operations 

dependent 960 540 60 90 

on Ns= 30 


Grand Total 1847 


1063 114 


12 96 


Blades , 

# Sq. 
Roots 


90 

91 
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today take an average of 1 to 6 microseconds per operation (large super 
computers excluded) the frame time would range between 23.6 and 142 
milliseconds, far more than our 4 millisecond goal for real-time RSRA 
simulation. And we haven' t included fuselage, wing, and tail aerodynamic 
forces in our RSRA model, which will add somewhat more time per frame. 

Assuming 50,000 normalized digital operations per pass as a rough 
overall estimate, including overhead, we conclude that a general purpose 
digital computer capable of 12. 5 million operations per second (50,000x 
250 hertz) is needed for the real-time RSRA simulation. The alternative 
is to use a specially configured digital computer with pipelining and some 
parallelism. Computers in this category with processing speeds in the 
correct range to solve the RSRA problem are, or soon will be, available 
at reasonable cost. Higher level compilers (e.g, , Fortran IV) are a long 
way off for such machines, however, so that assembly language programm- 
ing will be required, at least for some time, if analysis shows that floating- 
point versions of such machines are fast enough, then clearly that is the 
preferred alternate. With the scaling techniques used to develop the 
equations in Appendix B we are also convinced that a 16-bit fixed point 
digital mechanization can also do the job, if that approach becomes 


necessary. 



4. CONVENTIONAL. HYBRID SIMULATION 


4. 1 Basic Approach 

From the previous section it is clear that the speed capabilities of 
small to medium-size general-purpose digital computers are at least an order 
of magnitude below the requirements for a real-time RSRA simulation. One 
possible alternative would be to use conventional all- parallel analog 
mechanization for the portions of the simulation involving the rotor frequencies 
(3 0 hertz) and digital mechanization for the low speed (rigid airframe) 
simulation. Unfortunately, reference to Table 3. 1 shows that roughly 95 percent 
of the computational load involves the rotor frequencies (all the N or Ns 
dependent operations). With a conventional hybrid approach this implies a 
very sizeable analog requirement. 

4. 2 Hybrid Requirement 

Let us assume that a typical analog console includes as many as 
64 multipliers. In Table 3. 1 the total multiple /divide count was based on 
digital mechanization. If we go back to Table 2. 1 where the number of 
multiply /divides was also tabulated assuming analog mechanization, we calculate 
the following multiply /divide count for N = 6, s =5: 

# multiplier s /divider s = 80 x6 + 22 x30 = 1140 
Assuming 64 multipliers per analog console, this implies 18 analog consoles! 
Also, the 96 two -variable functions represent a formidable task if implemented 
with analog mechanization. To be sure, these could be generated using an 
available specially- configured digital computer at a frame rate sufficiently 
high (e. g. , ^ 500 hertz or above) to avoid degrading the hybrid simulation. 
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But if a specially-configured digital computer is used for the function generation, 
why not use it for the whole simulation, assuming the speed is adequate? 

Thus the conventional hybrid approach appears quite unattractive to us. 

There is one approach which might be considered. That alternative 
is to do all the computation except the actual integration on a specially- 
configured l6-bit digital computer with the integration performed on one 
analog console. This avoids the problems which would otherwise arise on the 
l6-bit computer if it were used for the numerical integration as well. I. e. , 
for the slow-varying state variables the incremental change in state from 
one frame to the next is too small compared with the state variable itself 
to permit reasonable sealing of the increment with a l6-bit computer. Having 
more than 16 bits (say 32) representing the integrands solves this problem 
but requires hardware development. Using analog integration is an immediate 
way around the problem. For the RSRA with 6 rotor blades there are 
6 + 12 = 18 degrees of freedom, or 36 integrators minimum. 



5. TIME-SHARED HYBRID SIMULATION 
5. 1 Basic Concept 

An alternative hybrid mechanization of the RSRA equations involves 
the use of an analog circuit to compute the aerodynamic forces and moments 
along each blade and, possibly, to solve the blade equations of motion. The 
single analog circuit is then time shared amongst the N blades. This approach 
has been the basis for hybrid mechanization of rotary-wing aircraft simulations 
in the past using general-purpose analog /hybrid computers, and is the concept 
used in at least one current special-purpose hybrid design for helicopter 
simulation. 

In this approach the aerodynamic forces on a given rotor blade are 
computed continuously on the analog computer as a function of distance along 
the blade, from hinge to tip. An analog integrator with fixed input provides 
a high-speed sweep voltage representing distance along the blade. The 
resulting aerodynamic force components at each point along the blade, as 
voltage outputs of the analog circuit, are integrated with appropriate weighting 
factors so that at the end of each sweep the total force and moment components 
acting on the blade are represented by the final integrator outputs. The 
integrators are then reset to zero and the sweep is repeated with the same 
analog circuit representing the next blade, etc. , until total force and moment 
components for all N blades have been computed. The process is then repeated 
cyclically, with the result that a descrete time-series of analog voltages 
representing force and moment components for each blade is generated. These 
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can then be used to solve the blade equations of motion, either with analog or 
digital integration. The result is a considerable reduction in required 
digital computer load for the overall simulation. By time -sharing a single 
analog circuit over all blades a substantial saving in analog equipment is 
achieved, 

5. 2 Equipment Requirements 

First of all, let us assume that the time -shared analog scheme described 

above is used only for computing overall blade force and moment components. 

Then the saving in digital operations applies to all the operation counts in 

Table 2. 1 involving s, the blade segment number. Thus the velocity- component 

calculation along the blades in Eqs. (8. 6)-(8. 8), all the aerodynamic force 

components in the equations of block 10, the blade moment computation in 

block 11, and the force summations in Eqs. (13. 6)-(13. 8), in short, all the 

equations in Table 2. 1 where the equipment count is multiplied by Ns, are 

included. Thus 32Ns multiplies, 18Ns adds, 2Ns trig functions, 3 Ns two- 

variable functions, and 3Ns square root digital operations per pass are saved 

by the time-shared analog mechanization. Clearly the latter requires 22 

analog multipliers (see the analog multiply /divide column in Table 2. 1), 

2 trig function generators, 3 two -variable function generators, and 3 square-root 

generators (analog squarers in amplifier feedback loops). This is a reasonably 

small complement of analog equipment, although the 2 -variable analog function 

generators might pose a problem. Actually, one of the three 2 -variable 

functions (the computation of OL a-nto) mechanized algebraically, so that 

-XRANb 
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only the 2-variable functions representing blade section lift and drag 
coefficients need be implemented. 

Assuming all this is done, reference to Table 3. 1 shows that for N = 6 
(6 blades) and s = 5 (5 blade segments in the equivalent digital mechanization) 
there is still considerable digital computation remaining, namely, the operations 
dependent on N and independent of N. These represent perhaps 40% of the 
original digital load for doing all the equations, so that the digital task is still 
formidable. Thus we have to perform more of the task on the analog, e. g. , 
integration of the blade equations of motion (Eqs. (9-5), (9-6), (12. 1), and 
(12. 2)), as well as a large number of coordinate conversions. In fact, if we 
are to make a major impact on the digital operations count, say reduce it by 
an order of magnitude, clearly most of the N dependent as well as all of the 
Ns dependent operations in Table 2, 1 must be performed on the time -shared 
analog circuit. In addition to the components required above for the Ns 
dependent equations (22 analog multipliers, etc. , ) reference to the summary 
at the end of Table 2. 1 shows that up to 80 analog multipliers and 8 trig 
function generators would be required. Although this is still much less 
analog equipment than would be needed for the conventional hybrid mechanization 
of Section 4, it would require a total of at least two analog consoles if general- 
purpose equipment were utilized. 

In the above time -shared hybrid mechanization the analog circuit 
inputs, coming from DAC^s (digital -to -analog converters), would be digitally- 
generated components of shaft-axis velocity, angular velocity, acceleration. 
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and angiilar acceleration. The analog circuit outputs, applied to A to D 
channels interfaced to the digital computer, would be a discrete time series 
of force and moment components along shaft axes due to each blade. After 
the time -shared analog circuit completes a duty cycle for a given blade, the 
force and moment contributions due to that blade are converted to digital 
form through the interface. The analog computation then proceeds through 
next blade duty cycle, etc. 

5 . 3 Dynamic Considerations 

From the previous section it is clear that the time- shared analog circuit 

is required to solve the blade equations of motion, i. e. , integrate Eqs. (9. 5), 

(9.6), (12. 1), and (12,2) in Appendix B. These integrations will proceed 

continuously over each individual blade duty cycle, but many of the terms on 

the right-hand side of the equations will be fixed over the duty cycle, e. g. , 

the blade moments m . and m . in Eqs. (9-5) and (9- 6). Thus the 

Pxj pzj 

mechanization for each degree of freedom (lagging and flapping) becomes in 

principle a second-order loop with two continuous analog integrators but with 

a sample -data feedback loop. The dynamics of such loops have been studied 

3 

extensively using z -transform theory . In particular, the damping ratio and 

frequency errors due to the finite frame time T of the feedback loop, plus the 

additional delay of T needed to calculate the quantity fed back (aerodynamic 

moment) turn out to be the following: 

3 

Damping ratio error - 

Fractional freq. error = - 0.864 (aT)^, a T << 1 (5-2) 

3 

Gilbert, op. cit. 
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Here T is the duty- cycle time for each blade computation (time to sweep 
from hinge to tip to calculate total blade moment components) and a is the 
frequency of the two -integrator loop. We will assume that this frequency 
turns out to be the order of the rotor spin frequency, i. e. , 30 rad/sec. 

Actually, the circuit should probably be mechanized with a separate 
pair of integrators for each degree of freedom for each blade, i. e. 4N 
integrators in all. Then the group of 4 integrators assigned to each blade 
will integrate only for 1/N of the overall duty cycle, and will be on hold for 
the remining time while the equations for the other N-1 blades are integrated. 

In effect this means that each analog integrator must be programmed to 
operate at N times real time, since it integrates over only 1/N of each total 
duty cycle (duration N times the individual duty cycle period, T). Thus the 
frequency a in Eqs. (5. 1) and (5.2) is effectively N times 30 or 3 ON rad/sec. 

If a damping ratio error of 0. 01 is deemed an acceptable upper limit, 
then Eq. (5.1) shows that aT = (4. 3)(.01) = 0.0133. This corresponds to 
75 frames per radian. For 6 blades (N = 6) a = 30(6) = 180 radians, so that the 
frame rate = 1/T = 75(180) = 13,500 frames per second. The corresponding 
duty cycle period T = 75 x 10 or 75 microseconds. This is probably 
unacceptedly short, since one would expect time -constants in the analog 
computing elements (especially the function generators) to be the order of 
one microsecond. Furthermore, mode control timing errors could be 1 micro- 
second. The result would be dynamic errors in the analog circuit of perhaps 
several percent. 
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Fortunately, the sampled-data feedback delay in each two -integrator 


loop is known (to first order it is 1. 5T) so that rate compensation for this 

delay can be implemented. The resulting errors in damping ratio and frequency 
2 

now depend on (aT) and acceptable dynamic accuracy (better than 0. 01 in 

3 

damping, one percent in frequency) can be obtained with aT = 0. 1. 

For a = 30(6) = 180 rad/sec, this corresponds to 1800 frames per second, 
or T = 556 x 10 ^ or 556 microseconds. This is a more reasonable duty cycle 
and should minimize the effect of analog component dynamic errors. 

To summarize, the time shared hybrid approach seems feasible for 
solving the RSRA equations, but will require at least two analog consoles. 
Naturally, it suffers from all the programming and reliability problems 
associated with such a hybrid approach. 


3 


Gilbert, op. 


ct. 
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6. SUMMARY RECOMMENDATIONS 


As a result of this study the following recommendations are made: 

1 . The real-time simulation of the RSRA equations should be 
implemented on a specially- configured digital computer using floating point 
if the speed is sufficient. If fixed point computation is necessary, it is 
recommended that equations similar to those in Appendix B of this report 
be used. 

2. It is recommended that NASA lead an effort to standardize the 
notation and sign conventions used in rotary- wing aircraft math models, and 
that this be as consistent as possible with conventional aircraft 

notation. 

3. It is recommended that NASA conduct an ongoing research program 
in rotary- wing aircraft modeling for real-time simulation. This effort should 
involve coordination of past and present efforts of the rotary -wing aircraft 
manufacturers, implementation of 2. above, systematic consideration of 
competing modeling approaches, determination of possible equation simplifications, 
investigation of different integration algorithms, determination of required 
integration frame rates, and correlation of results with flight tests. 
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APPENDIX A 


DEVELOPMENT OF HELICOPTER EQUATIONS OF MOTION 

AND ROTOR MODEL 


A. 1 IN TRODUCTION 

In this appendix we develop from first principles the helicopter 
equations of motion, including the rotor model. The purpose of this 
development is to verify the equations furnished by AUC, ^ understand 
certain simplifications made in those equations, and in some cases sug- 
gest alternate mechanization which may be more suitable for computer 
simulation. In Appendix B, the equations are rewritten in dimensionless 
form in order to facilitate the simulation if hybrid or fixed-point digital - 
computation is utilized. 

A. 1.0 NOTATION FOR VECTOR COMPON E NTS 

In our equations of motion we are working with components of 
vectors along a set of mutually orthogonal right-handed axes. The vec- 
tors represent velocities or accelerations of one axis system relative to 
another. For convenience we will designate the various axis systems 
with a single letter; thus the I-frame is an inertial reference system, 

B -frame is a set of body axes, etc. The individual axes will be desig- 
nated with a lower case subscript; thus x^, y^, are the x, y and z 
body axes. A vector velocity of one axis system relative to another 
usually will be designated with a double subscript; for example, is 

the angular velocity of the P-frame relative to the Q-frame. This vector 

^Mackie, D, Brian, and Thomas S. Alderete, ‘*A Real-Time Dual Pro- 
cessor Simulation of the Rotor System Research Aircraft, ** NASA TN- 
D-8328, January 1977. 
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can, of course, be projected geometrically into a third reference frame, 
say the R-frame. This third reference frame will be designated with 

a lower case superscript, where necessary for clarity. For example, 

- r r r 

the X , y , z components of will be written P , Q , R , 
r r r ^ pq pq pq pq 

respectively. 

We have used the convention that if the vector is relative to iner- 
tial space, the I-frame, then the second subscript may be omitted. 

Thus the R-frame components of ^2^^, the angular velocity of the P- 

r r r 

frame relative to the I-frame, may be written P , Q , R . If a single 

P P P 

subscript is used, it always is implied that the second subscript is ”i'*. 

In many cases the vector is projected onto the axis system given 

by the first subscript. When this is the case, the superscript may be 

omitted. Thus P , Q , R are the P-frame components of Q , and 
pq pq pq ^ pq 

P^, Q^, R^ are body axis components of 

A. 1.1 DEFINITION OF AXIS SYSTEMS USED 
A. 1. 1. 1 Inertial Axes, I- Frame 

For this study it is assumed that the earth is flat and non -rotating 
so that a frame translating with a constant velocity (i. e. , the horizontal 
wind velocity) with respect to the earth may be considered Newtonian. 

We designate this atmospheric frame as the I-frame. 

A. 1. 1. 2 Body Axes, B-Fram e 

These are conventional vehicle body axes. Origin of the B-frame 
is the vehicle center of mass. If there is a plane of symmetry (such as 
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most airplanes have) it is the plane. In straight and level flight 

the positive axis points approximately forward, the positive y^ axis 
points approximately to the pilot's right and the positive axis points 
approximately downward. The B -frame usually is defined by the vehicle 
manufacturer. 


A.. 1 . 1.3 Hub Axes, H-Frame 

Origin of the rotor-hub axes is displaced from the vehicle center 

of mass by the vector r , where 

JH 

V 

Hub axes are parallel to body axes, as shown in Figure 1.1. 


A. 1. 1. 4 Shaft Axes, S- Frame 

Origin of the S-frame is the same as hub axes. The S-frame is 
rotated from the hub frame through the fixed Euler angles i^ and i, . 

C7 (p 

The Euler angles result from positive rotation of i^ about the y^ axis 

followed by i, about the resulting x axis. 

cp s 


A. 1. 1.5 Rotating Shaft Axes, "Frame 

th 

Origin of the R^ -frame for the j blade is the same as the S-frame, 

namely the hub. As can be seen in Figure 1. 1, the x^^, y^^, 2 ^^ ^ R^- 

frame axes are rotated with respect to the x , y , z S-frame axes 

^ s s s 

by the azimuthal angle about the z axis. Thus the z . and z 

^ ® ^ j s rj s 

axes are coincident. The angular velocity of the R^ frame with respect 
to the S-frame is given by and is directed along the z^ (i.e., 
axis. Clearly There are separate R^ -frame axes for each 
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Figure 1. 1 Relative orientation of body axes (x , z^), hub axes 

(x^, y^, z^), shaft axes (x^, y^, z^), rotating shaft axes 
(x y z .) and blade “Span axes (x y z .) for the 

PJ PJ PJ 
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of N rotor blades, i. e. , j=l,2, 


Ao 1. 1.6 Blade Span Axes, Frame 

Origin of the P^-frame is at the blade -hinge point shown in 

Figure 1. 1. The hinge-point P^ is located a distance e along from 

the hub. In Figure 1. 1 there is indicated an axis system y^^, 

which is just the R. -frame axes, x V -i z translated from the hub 
3 rj ^rj rj 

to the hinge -point P^. The P^ -frame is rotated from the -frame 

through the blade lag -angle 5. (negative rotation about z* .) and the blade 

3 ^3 

flapping-angle j3. (negative rotation about x .). The y . axis of the P.- 
3 P3 P3 3 

frame is along the spanwise direction of the blade. Thiere are separate 
Pj -frame axes for each of N rotor blades. 


A. 2. HELICOPT ER EQUATIONS OF MOTION 
A. 2.0 SUMMARY 

A block diagram of the overall equations of motion in vector form 
is shown in Fig. 2. 1. It will be used as a reference for the listing of the 
scalar equations in the following sections. The latter are the basis for 
establishing the computational requirements for simulation. The blocks 
in Figure 2. 1 are numbered in the lower left or lower right corner. 

Block 1 represents the overall aircraft translational and rotational equa- 
tions of motion as determined by summing forces and moments about 
the aircraft, e.g. , Blocks 2 through 5 represent the equations for the 
accelerations and velocities of the rotor hub axes and rotating rotor 
shaft axes. Blocks 6 and 7 represent equations for the velocity and 
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Blod# Span Axit Acctl. Blodi Eqi. of Motion 



Figure 2. 1 Block Diagram of Helicopter Equations in Vector Form 














acceleration of the rotor blade -span axes. Block 8 shows the equations 
for the velocities of the individual blade segments. Block 9 represents 
the blade equations of motion, with two degrees of freedom (lagging 
and flapping) for each blade. Block 10 represents the equations for 
generating aerodynamic forces on the blade, while block 11 shows the 
equations for calculating the total external moment on the blade (except 
gravity). Block 12 represents the equations for relative velocity of 

blade span axes with respect to rotating shaft axes, which are used to 

. • 

establish equations for 6 and jS , where 6 and |3 are blade lag and flap 
angles, respectively. Finally, blocks 13, Ik, and 15 show the equa- 
tions for computing the overall vehicle force and moment due to the 
rotor forces » 

A. 2. 1 TRANSLATIONAL A ND RO TA TIONAL AIRCRAFT EQUATION S 
The translational equations for the aircraft are based on the vec» 
tor equation 

- G = F^/M (1. 1) 

Here A, . is the translational acceleration of the aircraft body axes 
bi 

(x^, y^, z^, centered at aircraft c.g. ) with respect to the inertial 
reference frame which is the atmosphere through which the aircraft 
flies. G is the gravity acceleration vector and F^ is the total of all 
other external forces acting on the aircraft. M is the aircraft mass 
(not including rotor blades). In scalar form, Eq. (1.1) becomes 


(1.3) 


A- - G = F, /M 
Dy y by 

'*■ “> 

Thus the external force components, ^by' ^bz ^sed to 

provide the three acceleration-less -gravity output components of block 
1 in Figure 1.1. We note that 

*bi = ^ ^ Vb ^ '‘-5) 

Here is the body axis translational velocity vector with respect to 

the inertial reference frame and ^2^^ is the body axis angular velocity 

vector with respect to the same frame. ""V. . is the time-rate-of- 
^ dt bi'i 

change of V . as viewed from the inertial (I) frame, and ~ V, .)_ is the 
bi at bi b 

time -rate -of-change of as viewed by an observer in the body axis 
(B) -frame. Thus 


^bi " Jk + 


b b 


b ^b 


b b 


( 1 . 6 ) 


where i^, and k^ are the unit vectors, respectively, along the x^, 

y, , body axes. Hence U. , V, , W, are the components of aircraft 
b b b b b 


velocity along the three body axes, respectively. Also 


= Pb\ + °bib + ^\ <‘-"i 

where P^, Q^, are, respectively, the components of aircraft angular 
velocity along the three body axes. Thus represents conventional 
aircraft roll rate, the pitch rate, and the yaw rate. 

From Eq. (1.6) we obtain 

^‘>b - “b‘b*^bi'b^*b\ 
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since it follows that for an observer in the body-axis frame, the body 
axis unit vectors i^, do not change. We also note that the 

gravity accelerations in Eqs. (1. 2) through (1. k) are given by 


G 

X 

= - cm 9^ 


(1.9) 

G 

y 

= g COS 0^ sin 

^o b 


(1. 10) 

G 

z 

= g cos 0- COS 

b 

‘t’b 

(1. 11) 


where 0^ and <(>^ are the conventional aircraft pitch and bank angle, 
respectively, and where g^ is the sea level gravity acceleration. 
Substituting Eqs. (1.6), (1.7) and (1.8) into (1.5), and the scalar com- 
ponents of (1. 5) into (1.2) through (1.^), we obtain the following three 
equations using (1.9) through (1. 11): 

' - "b°b * ''b'H, - So «b <>• 

\ ' •''b'H, '^b‘“b * cos 9^ sib (1.13) 

*b = - ''b^’b “b°b * So <=■=■' ^ +b * ''bz/“ 

These translational equations of motion are integrated to obtain the body- 
axis velocity components U^, V^, and W^. 

The rotational equations of motion are obtained from the following 
equation: 

^'i = “b 

where is the external moment acting on the aircraft c. g. and is 
the angular momentum of the aircraft about the c.g. is in turn 

given by 
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( 1 . 16 ) 


Here [l] is the inertia matrix for the aircraft, not including the inertia 
of the rotor blades. This is because the angular momentum associated 
with the blades is considered separately in the blade equations of motion, 
and the resulting blade reaction moments on the hub are part of the total 
external moment acting on the aircraft. Also in Eq, (1. 16) we have 
not included the angular momentum of the turbine power-plant assembly, 
although this can easily be done. 

If we neglect products of inertia in the inertia matrix [l], then Eq, 
(1. 16) can be written as 



where I , I and I are the aircraft moments of inertia about the body 
XX yy zz ^ 

axes, x^, y^, z^, respectively and i^, and are the unit vectors 
along the respective body axes. Differentiating Eq. (1, 16) with respect 
to time we have 


* “bi * 

n Eq. (1. 18), ^ ^bi^b ' the time rate of change of body-axis angu- 

ar velocity as viewed by an observer in the body axes, is given by 


dt ^bi^b ^b ^b ^b ^b ^b ^b 


(1. 19) 


ince we have neglected the products of inertia. 




( 1 . 20 ) 


34 



and hence 


’ «zz-V°bVb ^ 

^ ' 1-2 

Substituting Eqs. (1.19) and (1 . 21 ) into Eq. (1.18), and (1.18) 
into (1.15), we obtain the following thress scalar equations of motion 
of the aircraft: 


I -I 
zz yy 

I 

Q £b + 

(1.22) 

XX 

XX 


I - 1 

XX zz 

I 

■R,P 

I 

(1.23) 

yy 

yy 


I -I 

M, 


yy XX 
I 

■ P Q, + — ~ 
b^b I 

(1.2lf) 

zz 

zz 



^y' ^bz 

represent the exter 


In Eqs. (1.22) through (1 
nal moments acting along the respective body axes. 

lii case the products of inertia are not negligible, it is straight- 
forward to add the requisite terms to Eqs. (1.22) through (1.24). These 
equations are integrated to obtain the body-axis angular velocity com- 
ponents, P^, Q^, which in turn are used to compute the time-rates - 

of-change of the body axis Euler angles (heading), 9 ^ (pitch) and cj)^ 
(bank angle). Thus we have the familiar formulas 



= (R^ cos cf)^ + sin (j)^)/ cos 6 ^ 

(1.25) 

"b 

= cos 4>b - R-b ‘t’b 

(1. 26) 


= P, +4!, sin 0, 
b b b 

(1.27) 
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These are integrated in turn to compute the Euler angles 0^, and 




If our simulation requires computation of the aircraft position, 

it is necessary to calculate aircraft velocity components along earth 

— 0 

axes, e. g. , north, east, and vertical. Denoting as the aircraft 

velocity vector expressed in components along the earth (e) axes and 

as the aircraft velocity vector expressed in components along 
the body axes, we write 

where [T^^] is the matrix used to transform body-axis components of a 
vector to earth-axis components. It is straightforward to show that 


[T 


be 


] = 


cos G 
cos G 


, sin^- 
b b 

- s in 0 , 
b 


-coscb^ sin^, +sincb, sin0, cos^ 
^b b b b b 

cosc^^cos ^+s incjj^s in^^sin"^^ 


s iiKb^ cos 0- 
^b b 


sin6,sin^- +cos<|), sin0_ cos^, 
b b b b b 

•»sinct)^cos^+co3 4>^sin0^sin^ 


cps coscb, 
b b 


(1.29) 


• • • . 

Letting S , S , and S (i.e., -n ) represent the velocity components 
X y 2 

north, east, and downward (h is altitude), we have from Eqs. (1.28) 
and (1.29) 


S 

X 


U, cos 9, cos^, 
b b b 


+ V, (-cos A sin + sincj), sin0, cos 
b b b b b 




+ W^(sin(j)^sin^^ + cos<j)^sin0^cos^) (1.30) 

S = U, cos0, sin^, + V, (cos6, cos^, + sincj), sinO, sin^, ) 
y bbob bb bbD 

+ W. (- sind), cos>fe + cos(|>_ sin0, sin^_ ) (1. 31) 

b ^b b b b b 
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S = - li = - U, sinO, + V- sin<(), COS0, + W, cosO. cos<b- (1. 32) 

z DbDbbbbb 

We note in passing that the inverse of the matrix is its transpose. 

Thus 


where 





(1.33) 


and is given in Eq. (1*29). 

In this section we have assumed as an inertial reference frame 

the atmosphere through which the aircraft is moving. If there are 

steady wind components they are easily added to Eqs. (1. 30) through 

(1. 32) to take care of wind drift. 

In summary, we have in this section indicated the derivation of 

the six-degree - of - freedom equations of motion of the rigid aircraft. 

Three of the state variables are the translational velocity components 

U^, and along the body axes. Equations (1. 12), (1. 13) and (1. 14) 

are integrated to obtain these variables. Three more state variables 

are the rotational velocity components Q^, along the body axes. 

Equations (1.22), (1.23) and (1.24) are integrated to obtain these 

variables. The Euler angles 9 ^, constitute another three state 

variables, and Eqs. (1. 25), (1.26) and (1.27) are integrated to compute 

them. Finally, distance north (S ), distance east (S ) and altitutde 

X y 

(h = - S ) are the remaining state variables. Equations (1. 32), (1. 33) 
z 

and (1.34) are integrated to compute them. 

Although the equations derived in this section are well known and 
widely used, it was felt worthwhile to include the derivation because the 
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methodology carries through to the subsequent derivation of rotor blade - 
equations of motion. 


A. 2.2 ROTATING SHAFT AXIS ANGULAR VELOCITY 

Although the rotor hub is located at a position remote from the air- 
craft c. g. , reference to Figure 1. 1 shows that the hub axes y^, 
are aligned with the body axes x^, y^, z^. Thus the angular velocity 
of the two axis systems is identical, and hence we write 

Next we determine formulas for the components of along the shaft 

axes X , y , z . From Figure 1. 1 we see that the shaft axes differ from 
s s s 

the hub axes by the pitch angle i and roll angle i . Denoting the trans- 

C7 tp 

formation matrix from hub to shaft axes as [T, ], we have 

ns 




( 2 . 2 ) 


“S 


where indicates that is expressed in terms of shaft-axis com- 

— h “ 

ponents and indicates that is expressed in terms of hub-axis 
components . 

The matrix is obtained easily by analogy with conventional 

Euler angles, for which we have already written the matrix in 

Eq, (1,29) for the transformation of vector components from body to 
earth axes. In Eq. (2. 2) is equivalent to a matrix for transform- 

ing vector components from earth to body axes where 6 ^ = i^, 

cj)_ = i., and = 0. Thus we set 0, = i^, c|), = i, , and = 0 in 
bcp b bC7b<p b 

T 

[T^^] = to obtain Making these substitutions in the trans- 
pose of i^ Eq. (1.29), we obtain 
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cos 1 , 


- sin 1 . 


sini^sini^ cos sini^cosi^ 


cosi.sini^ - sini, 
<p y <p 


cos 1 , cos 

9 y 


(2.3) 


Since and are fixed angles, the angular velocity of the shaft axes 
Xs* y ^9 Zg is equal to the angular velocity of the hub axes x^, y^, z^. 
Hence, we can write 

“=l = <2. 4) 

Thus, Eqs. (2.2) and (2.3) can be used to compute the components of 

shaft-axis angular rate along the shaft axes- 

Reference to Figure 1. 1 shows that the rotating shaft axes x 
th 

y z .for the j rotor blade differ from the shaft axes x , y , z 
rj r j s s s 

by the azimuthal angle in the x^, y^ or x^^. y^^. plane, where 


= -a 

J s 


(2.5) 


Note that and hence the rotating shaft axes are different for each of 
the N blades (j=l, 2, . . . , N). To obtain the angular velocity of the 


rotating shaft axes for each blade, we note that 


Q . = d , + n k ( 2 . 6 ) 

ri SI s s 

where is the unit vector along the z^ axis and is assumed to be the 
velocity of the rotating shaft axes with respect to the shaft axes (see 
block 2 in Fig. 2^1). Note that the angular velocity of the rotating shaft 
axes is the same for all blades. This is why we denote the angular ve- 
locity as (angular velocity of R. -frame with respect to the inertial 

th 

frame) without an added subscript j for the j blade. 
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Using Eqs. (Z. 2), (2.3), (2.4) and (2.6), we obtain the following 

s s s 

three scalar equations for the components and of rotating 

shaft axis angular velocity along shaft axes: 


P = P, cos - R, sin i^ 
r b 0 b 0 


Q = P, sin i. sin i^ + Q, cos i, + R, sin i, cos i^ 
r bc|)0b 0 


(2.7) 

( 2 . 8 ) 


R = Pi_ cos i, sin i^ - Q, sin i , + R_ cos i, cos i^ + Q 
r b <p0b(pb 9 8 i 


(2.9) 


s s s 

Note that P j Q , and R are the same for each blade. For this reason 
r r r 

th 

we don^t add the subscript j for the j blade. 

Next we determine the components of along the rotating shaft 

axes x^, y^, z^. Again, reference to Figure 1. 1 shows that the rotating 

th 

shaft axes x y ., z , for the j blade differ from the shaft axes x , 
rj ^rj rj s 

y , z by the equivalent Euler angles 0_ = 0, 4), = 0, ^ = tt/ Z - 

ss b b b j 

T 

Making these substitutions into = ^^be^ ' where is given 

in Equation (1.29), and noting that cos(;r/2 - = sin 

sin {tt/Z - ^) = cos we obtain the following transformation matrix: 


[T ] = 

sr/ 

J 


sin cos 0 
J J 


- cos sin 

3 3 


( 2 . 10 ) 


Thus 




ri. 

J 


[T ]S2®. 
sr. ri 
J 


( 2 . 11 ) 
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and from (2. 10) and {2, 11) we have 


r s s 

P , = P sin + Q cos 

rj r j r ] 


Q^. = - P^ cos sin 

rj r j r j 

R . = R 
rj r 


( 2 . 12 ) 

( 2 . 13) 
(2. 14) 


as the components of rotating shaft axis angular velocity along the R. 


J 


.th 


frame axes for the j blade. Since there is a rotating shaft axis sys- 
tem for each blade, Eqs. (2. 12) through (2. 14) are repeated N times 
for N blades, each blade having a different (j = l, 2, . . . ,N). 

Note that the N blades will be equally spaced in azimuth with a 

separation of 2 7 t/N between each blade. Thus if is the azimuthal 

th 

angle for the first blade, then the angle >3^ for the j blade is given by 


'Jr. = + 2jr(j-l)/N, 3=1.2,..., 


N 


(2.15) 


It follows that 


sin = sin[^ + ^ 

3 1 N 


and 


sin cos 




( 2 . 16 ) 


cos ^ = cos[^^ + = cos cos[- — - sin sin [ 

(2. 17) 

Thus once sin and cos are computed for the first blade, sin ^ 
and cos for the blade can be computed using Eqs. (2. 16) and 
(2. 17). Each equation only requires two multiplies and an add, which 
may save considerable time over computing sin and cos ^ directly. 
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Finally, we develop the equations for obtaining the components of 


rotating shaft axis angular velocity along the blade span axes for the j 


.th 


blade. In Figure 1. 1 we see that the blade span axes, x y z . are 
® ^ PJ PJ PJ 

displaced from the rotating shaft axes x ., y and z . by the yaw angle 

b. (called the blade lag angle) and the roll angle (called the blade flap 

angle). Again, using the Euler angle analogy, 0^ = 0, = - p^, and 

- 5 ., Making these substitutions into [T^] = [T, ] , where [T- ] 

b j eb be be 

is given in Eq. (1. 29), we obtain the following formula for the transfer- 

th 

mation matrix from rotating shaft axes to the j blade span axes: 


[T 




since 


we have 


cos 5 


J 


sin 5 ^. 


sm 5 . cos B. cos 5 . cos 5 . - sin q. 

J ^3 3 ^3 ^3 


sin 5 . sin 6. cos 5 . sin b. cos S- 

3 3 ^3 ^3 J 


= [T 

ri. rp ri. 

3 Pj J 


(2. 18) 


(2. 19) 


p^. = 

M 

cos 

5. 

- 

sin 5 . 




( 2 . 

20 ) 

rj 

^3 


3 

^3 

3 






Q^. = 

= P^. 

sin 

5 . 

cos 6- 

X 

+ Q . cos 

5; 

^r 

cos p - R 

sin R 



rj 

^3 


3 

^3 

^3 

3 

^3 ^3 

J 












(2. 

21) 

rP. = 

= 

sin 

6. 

sin p. 

r 

+ Q . cos 

5. 

r 

sin B. + R . cos B. 

(2. 

22) 

rj 



3 

^3 

^3 

3 

^3 ^3 

3 




Actually, it turns out that in Eq. (2. 21) is not needed in the mechan- 
ization. 

Thus block 2 in Figure 2. 1, which involves the computation of 
rotating shaft angular velocity components, requires the use of Eqs. 
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(2.7), (2. 8) and (2. 9); Eqs. (2.12), (2. 13) and (2. 14) times N for N 
blades; and Eqs. (2.20) and (2.22) times N for N blades. 


A. 2.3 ROTATING SHAFT AXIS ANGULAR ACCELERATION 

r_ 

Next we consider the equations for the angular acceleration 

of the rotating shaft axes x , y , z*. Consider firsts the angular 

r r r si 

acceleration of the shaft axes x , y , z • Since the alignment of these 

s s s ® 

axes is fixed with respect to the body axes (the angles i^ and i^ are con- 




stant in Figure 1,1), clearly 




%i 


(3. 1) 


where is the angular acceleration of the body axes. In section A. 2. 1 
we noted from Eq. (1.7) that the body-axis angular velocity, is by 

definition given by 


^i = ^b\ 

Differentiating we have 


dt ^bi\ dt ^bi^b ^i ^ ^Di ^ dt ^i^b 


(3.2) 


(3.3) 


since x = 0. Equation (3.3) shows that the time - rate -of-change 
of as viewed from the inertial frame is the same as the time rate of 
change as viewed by an observer in the body axis frame. But from Eq. 
(3.2) this is simply 




(3.4) 


where is the time rate of change of as viewed from the inertial 
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frame. Thus the body axis components of the output from block 1 

• • • 

in Figure 2. 1 are simply and respectively, as given in 

Eqs. (1. 22), (1, 23) and (1.24). 

The shaft axis components of can be obtained from body axis 
components by using the transformation matrix [T^g] = ^'^bs^ given 
in Eq. (2.3). Thus 


n®. = [T, ]n^. 

si hs SI 


(3.5) 


and since from Eq. (3, 1), ^2 . = we have 

^ si bi 


p® 

s 

II 

cos i^ - R^ sin i^ 

6^ 

s 


sin i^ sin i^ + 

R® 

s 


cos i^ sin i^ - Q, 

q? C7 D 


(3.6) 

(3.7) 


as the formulas for shaft-axis angular acceleration components along 

shaft axes. Next we compute the components of along the rotating 

th ^ 

shaft axes of the j blade using the formula 


= [T 

SI. sr. SI 

J J 


(3.9) 


where [T ] is given in Eq. (2. 10), Section A. 2. 2. Thus 

P^. = P® sin + Q® cos 'S', 

sj s J s J 


(3. 10) 


• r • s • s . 

Q . = - P cos + Q sin 

SJ s J s J 

;^s 

R . = R 

SJ s 


(3.11) 
(3. 12) 


We have already seen in Eq, (2. 6) of Section A. 2. 2 that ^2^^, the angular 


velocity of the rotating shaft axes, is given by 
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(3.13) 


^2 . = k 

ri SI s s 

Differentiating we obtain 


^2. = = -^^2.) +^2.xn, = -t:^2.) (3. 14) 

ri dt ri 1 dt ri r ri ri dt ri r 


since ^2 . x J2 . = 0. Thus from (3, 13) and (3. 14) 
ri ri 

j 

^ • = “jr(^ • + ^ ^ ) 

ri dt SI s r r 


(3. 15) 


where “s’ though there are N rotating shaft axes (one for each 

of N blades) until now we have not designated either ^2 . or S2 . with an 

® ri ri 

th 

additional subscript j for the j blade because each rotating blade-axis 
differs from the other rotating blade axes by a fixed azimuthal angle. 


Hence ^2 . and S2 . are the same for all blades. But as soon as we write 
ri ri 

the formulas for the components of t2^^ along the rotating shaft axes for 
th 

the j blade, we must use the additional subscript j. Thus in Eq. (3. 15) 
for the blade 


— r r - r - r - 

f2. = P.i+Q.j +R.k 

si^ sjr sj r sj r 

and from Eqs. (3. 10) through (3. 16) we have 


(3. 16) 


. P . = P sin ^. + Q cos 

rj s J s J 

• r • s • s 

Q . = - P cos ^. + Q sin ijr. 
rj s 3 s J 

• r • s 

R . = R + n 
rj S3 


(3. 17) 
(3. 18) 
(3.19) 


To summarize, the components of rotating shaft axis angular accelera- 

th 

tion along rotating shaft axes for the j blade are first obtained from 
Eqs. (3.6), (3.7), and (3. 8) followed by Eqs. (3,17), (3, 18) and (3. 19). 


Actually, in Eq. (3. 18) is not needed. Equations (3. 17), (3, 18) 
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and (3. 19) must be repeated N times for the N blades. 


A. 2.4 HUB AXIS VELOCITY 

The computation of hub -axis velocity, is indicated in block 4 


of Figure 2. 1. It is given by the formula 


V = V + V 
hi bi hb 


(4.1) 


where is the velocity of the body axes (i. e, , the aircraft c. g. ) and 

V is the velocity of the hub with respect to the aircraft c. g. (body- 
nb 


axis origin). is given by 




(4.2) 


where is the body-axis angular velocity and r^^ is the position of the 
hub with respect to the aircraft c.g. From Figure 1. 1 we see that 

where i^, j and k^ are the respective unit vectors along the body 
axes. Noting that ^ ^b b ^ ^ ^ have from Eqs. (4. 2) 

and (4.3) the following formulas for the body-axis (and hence hub-axis) 
components of 


U 




hb 

^hb " ^ 


^b^H 


^hb ^b ■ ^b 


(4.4) 
(4. 5) 
(4.6) 


To get the components of along shaft axes x^, y^, z^, we need 
to use the transformation matrix given in Eq. (2. 3). Thus 
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(4.7) 


= [T, ]vj. = [T, ](V^. + ) 

hi hs hi hs hi hb 


Finally, we note that 


V . = V, 

ri hi 


(4.8) 


i.e., the translational velocity of the rotating shaft axes is the same as 
the hub axes, since they have a common origin. From Eqs. (2. 3), (4. 7) 
and (4. 8), we obtain 






+ (W, + W, , )sin i. cos 

b nb q> y 


(4.9) 


(4.10) 




<“b ‘4. ‘e - <’^b * ''to’ S 


-t (W, + W, , )cos i cos 

b hb 9 o 


(4. 11) 


Fqs. (4.4), (4. 5) and (4.6), along with Eqs. (4.9), (4.10), and (4.11), 
are used to compute the velocity outputs of block 4 in Figure 2. 1. As 
before, we note that i^ and i, in Eqs. (4. 9), (4. 10) and (4. 11) are con- 

U cp 

stants, so that the coefficients in the equations can be precomputed as 
fixed parameters. 

A. 2. 5 HUB A XIS ACCELERATION 

The computation of hub -axis acceleration, A is indicated in 
block 5 in Figure 2. 1. It is given by the formula 




(5. 1) 
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Biit from Eq. (4. 1) V_ . = V . + V, , . Substituting (4. 1) into (5. 1), and, 

ni no 

noting that 5 “^i' have 

or 

K- = + "TT ^ ■^i,K (5.2) 

nx bi dt hb b bi hb 

But from Eq. (4, 2), ^ since r^ is fixed in the body- 

axis frame. 


dt ^hb^b ~ dt ^^bi^b ^ ~ ^ 

where 

" d d d 

bi dt bi 1 dt bi b bi di dt bi b 

Substituting Eq. (5. 3) into (5. 2) and subtracting the gravity acceleration 

vector G, we obtain 


° ° * 'h * \b '5-^> 

which is just the equation shown in block 5 of Figure 2. 1. Writing Eq. 
(5. 4) as three scalar equations, we obtain 


^z-°* = - °x + “b - 'S, * °b “h - ''h 

V- °y- ' ^b * ^b ' K "h 


As before we can transform these body-axis (and hence hub-axis) 
components of - G to shaft-axis components using the transforma- 
tion matrix [T, ] in Eq. (2. 3). Thus 

ns 
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(5.8) 


- 5 ”) 

We also note that the hub axes in Figure 1, 1 have a common 

origin with the rotating shaft axes x^, y^, z^. Hence 


From Eqs. (2.3), 
scalar equations: 


A . - G = 


^i “ 


(5.9) 


(5, 8), and (5. 9), we can write the following three 


K. °:=> ■=“* 'e - »“ ‘9 '=• “» 

s s 

A - G = (A, - G ) sin i, sin i^ + (A, - G ) cos i, 

ry y hx x <p 0 hy y <j> 

+ (A, - G )sin i. cos i,, (5. 11) 

hz z 9 C7 

s s 

A - G = (A, - G )cos sin i^ - (A, - G )sin i, 

rzz hxxcf>0hyy<t) 

+ (A- - G ) cos i, cos i^ (5.12) 

nz z 9 

As before, the coefficients can be precomputed as fixed parameters 

since i, and i^ are constants. 

9 tf 

Equations (5. 5) through (5. 7) and (5. 10) through (5. 12) are used to 
compute the hub-axis and hence rotating-shaft-axis acceleration com- 
ponents along the shaft axes. 


A. 2. 6 BLADE SPAN AXIS VELOCITY 

In Section A. 2.4 we determined the formulas for calculating the 

rotating-shaft-axis velocity components and alone the shaft 

r r r ® 

th 

axes. To obtain the components along the rotating shaft axes for the j 
blade we use the formula 



[T 


s r. 


J 



( 6 . 1 ) 
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where the transformation matrix [T ] is given in Eq. (2. 10), Thus 

s r , 

J 

from Eqso (2. 10) and (6. 1) we obtain 


U^. = U® sin + V® cos 

rj r J r j 


( 6 . 2 ) 


r s s 

V . = - U cos iSf. + V sin 

^3 r j r j 

W^. = W® 
rj r 


(6.3) 

(6.4) 


We next determine the formula for the translational velocity, , of 
th 

the j blade hinge point in Figure 1. 1. This is given by 


^ . = V . + V 

pi. ri. pr. 

J 3 ^3 


(6.5) 


where V is the velocity of the hinge point P. relative to the hub (i.e., 

p^j J 

the origin of the rotating shaft axes). Thus 


V = n . X e. 
pr. ri 3 


( 6 . 6 ) 


Here U is the angular velocity of the rotating shaft axes (we recall that 

. is the same for all blades, i. e. , is independent of j) and e. is the 
ri 3 

position of the hinge point P^ with respect to the hub. Thus 


e. = ei . 
3 rj 


(6.7) 


where is the unit vector along the y . axis in Figure 1. 1. Noting that 

n . = P^.T .+Q^.j .+R^.k . andV. = U^.i . + V^.j . + W^.k 
ri rj rj rj-'rj rj rj ri^ rj rj rj rj rj rj 

we have from Eqs. (6. 5), (6.6) and (6. 7) 


r r r 

U . = U . - R .e 

PJ rj rj 

V^ = v". 

PJ rj 

W^. = W^. + P^.e 

PJ rj rj 


( 6 . 8 ) 

(6.9) 

(6. 10) 
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Thus Eqs. (6.2) through (6.4) together with Eqs. (6.8) through (6. 10) 

are used to compute the rotating shaft axis components of the hinge ve- 
th 

locity for the j blade. All equations must be solved N times for N 
blades and provide the output of block 6 in Figure 2. 1. 


A. 2.7 BLADE-SPAN AXIS ACCELERATION 

The blade span axis acceleration A . is the acceleration of the 

pi. 

th 

hinge point P. for the j blade with respect to the inertial reference 
(I-frame) . It is computed using the same approach as in Section A. 2. 5 
for the hub acceleration. Thus by direct analogy with Eq. (5. 4) 

A. - G=A.-G+^2.xe.+^2.x^2,xe. (7.1) 

pij ri ri j ri ri j 

Here A^. = A^. (the rotating -shaft axes and hub axes have a common 

origin), is the angular velocity of the rotating shaft axes, and e^ is 

the position of hinge-point P. with respect to the hub. . x e. = V 
_ ^ ^1 J ptij 

= V , the velocity of hinge-point P. with respect to the hub. 

prj J 

Let us write Eq. (7. 1) in terms of scalar components along the 

th 

rotating shaft axes, x y ., z . for the j blade. To do this we first 

rj ^rj rj 

need to write the formulas for R-frame components of A^^ - G in terms 
of S -frame components using the transformation matrix [T ] given in 
Eq. (2. 10). Thus 


A^ . 
rxj 

- G^. 
xj 

= (A® - g' 

rx 

A^ . 
ryj 

1 

0 

= - (A® - 

rx 

A^ . 
rzj 

-g". 

zj 

= A® - G® 
rx z 


s s 

- G ) cos 

ry y J 


G^)cos 'ir. + (A^ - G^)sin 

X j ry y J 


(7.2) 

(7.3) 

(7.4) 
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where (A® - G®). (A® - G®), and (A® - G®) are given in Eqs. (5. 10), 

rx X ry y rz z 

(5. 11) and (5. 12). 

The R-frame components of U are given in Eqs. (3. 17), (3.18) 
and (3. 19); e^ is given in Eq. (6. 7); the R-frame components of 
are given in Eqs. (2. 12), (2. 13) and (2. 14). These equations provide 
all the R-frame components necessary to express Eq. (7. 1) as three 


scalar equations 

along the 

resp 

ective 

rotating shaft axes 

y z .. 

rj rj 

Thus we obtain 
A^ . 

- G^. = 

A^ 

- G^. 


• r 

- 

(7.5) 

pxj 




rj rj 



A^ 

- G^. 

A^ 

- G^. 

- [(R^)^ + 


(7.6) 

pyj 

yj 

""yj 

yj 




A^ , 

- G^. = 

A^ 

- G^. 

+ (o’". R^". 

+ P^)e 

(7.7) 

P2J 

zj 


zj 

rj rj 

r 



Equations (7,2) through (7. 7) must of course be repeated N 
times for each of N blades (j;=l, 2, . . . , N). 

A. 2. 8 BLADE SEGMENT VELOCIT Y 

In order to compute the aerodynamic forces acting on each of 
the blade segments, it is necessary to compute the velocity vector 
^jk center of pressure of the k*^^ segment of the blade. 

This is simply given by the vector formula 
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( 8 . 1 ) 


V., = V . +S2 . xp., 

Jk Pij Pij 

where V . is the velocity of the blade hinge at point P. in Figure 1.1. 
Pij 3 

^ is the angular velocity of the j blade -span axes, and is the 

J th 

position of the k segment center of pressure with respect to the hinge 


point Pj- Thus 




( 8 . 2 ) 


where y is the distance from P. to the segment center of pressure and 
Z K j 

is the unit vector along the y^^ blade span axis in Figure 1.1. We 



denote the blade -span -axis angular velocity the j blade by 


.. = P ,i . + Q J . + R 

P’3 PJ PJ PJ PJ PJ Pj 


(8.3) 


The components of blade -span -axis velocity V . have already been given 
in Eqs. ( 6 . 8 ), (6.9) and ( 6 , 10). To resolve them into components along 


the blade span axes x ., y ., z . requires the transformation matrix 

PJ PJ P3 

[T ] given in Eq. (2. 18). Thus 

''■j 


= [T 

pi. rp. pi. 

J J J 


(8. 4) 


or from Eq. (2. 18) 


p r r 

U . = U . cos 8 . - V . sin 5 . 


PJ 


vP. 

PJ 


PJ J PJ J 

r 


U sin 5 cos R. + V . cos 5 . cos R. - W . sin r. 
PJ J ^J PJ J ^J PJ ^^J 


p r r r 

W . = U . sin 8 . sin e. + V . cos 8 - sin 6 . + W . cos 6 . 

PJ PJ J '^J PJ J J PJ J 


(8.5) 

( 8 . 6 ) 

(8.7) 


Substituting Fqs. (8.2) eind (8.3) into Eq. ( 8 . 1 ). we obtain the following 
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the tangential 


three scalar equations for k' ^P;jk 

radial, and perpendicular components of the segment velocity 
along the respective blade span axes: 

u^., = U^. - R . y^, 

Tj.k pj pj ^2k 


= V^. 
P)k PJ 


w ■ + P 70, 

PJk PJ PJ 2k 



( 8 . 8 ) 
(8.9) 
( 8 . 10 ) 


Equations (8* 5) through (8. 10) are mechanized in block 8 in Figure 2. 1. 
Equations (8. 8), (8. 9) and (8. 10) must be repeated Ns times for N 

blades (j = l, 2, - . . , N) and s segments per blade (k = l, 2, . . . , s). 

A. 2,9 ROTATIONAL EQUATIONS OF MOTION FOR THE BLADE 


The rotational equations of motion of the blade are obtained by 

summing moments acting on the blade about the hinge point P. in Figure 
— th 

1.1. If we let H . be the j blade angular momentum (with respect to 


the inertial frame) about the point P., M^. be the total external moment 

J Bj 

about P acting on the blade, and A . be the acceleration (with respect to 
j Pi- 

2 

the inertial frame) of the point P^, then we can write 


H . ). = M^. - p , X m A . 

dt pu 1 Bj cj pu 


(9. 1) 


Here p is the position of the blade c. g. with respect to the hinge point 

P. and m is the blade mass. Thus - p . x m A . is the additional iner- 
J CJ Plj 

tial moment due to the acceleration A . of the point P., 


2 

Greenwood, D. T. , "Principles of Dynamics, Prentice-Hall, Inc., 
1965, p. 146. 
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Now H . = . » where [l^] is the blade inertia, matrix and 

pi B pi B 

— J J th 

Q . is the cingular velocity of the j blade. We will assume that 

0 0 




0 


0 

0 




(9.2) 


i.e., the blade has equal moments of inertia, L , along the x . and z . 

^ PJ PJ 

blade -span axes, and negligible inertia along the blade -span direction. 


the y . axis. From Eq. (9.2), then 
P3 


pi. D 

J 


P .i . + L R . k . 
P] PJ D pj pj 


(9.3) 


where P . and R . are the blade -span -axis (and hence blade) angular 
PJ PJ 

velocity components along x . and z respectively. Note that Q 

PJ PJ PJ 

the blade -span- axis angular velocity component y^^, is not equal to the 

blade angular velocity component along that axis, since in general the 

blade can be twisted (pitched) with respect to the blade -span-axes. 

But since the blade inertia about the y^^ axis is neglected, we don't 

need to consider the y . angular velocity component in representing the 

PJ 

blade angular momentum H . . 

Next we note that 


= +S2. xH. 

dt pi^ 1 dt pi^ p pu pi^ 


(9.4) 


where 


^2. =P.i.+Q.j.+R.k. 

Pij PJ PJ PJ PJ PJ PJ 

From Eq. (9.3) we see that H . ) is simply 

dt pi. p 

J 


(9.5) 
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(9.6) 


d t; 


^H.) = LP.i.+L 

dt pi, p ^ pj pj pj p3 


R . k . 


Substituting Eqs. (9.3), (9.5) and (9.6) into ((.4), we obtain 


-^H 
dt pi 


. ). = L (P . + Q .R ,)i . + L (R . - P .Q .)k . 
1 . 1 b PJ PJ PJ PJ b pj pj p3 pj 


3 


(9.7) 


Next we let M . represent the external moment exclusive of grav- 
PJ 

1 .th , , _ 

ity acting on the j blade, e-g. , 


M_. = M . + p . X m G 

B1 PJ cj 


(9.8) 


Substituting Eq. (9.8) into Eq. (9. 1), we obtain 


H.). = M.-p.xm(A. -G). (9.9) 

dt pi. 1 PJ pu 

But we have already developed in Section A. 2, 7 the equations for 
(A . - G), except that in Eqs* (7.5), (7.6) and (7.7) the components 

are along the rotating shaft axes. To convert them to components along 
the blade -span axes the transformation matrix [T ] in Eq. (2. 18) is 
used. Thus 


iP - 


(A . - G ,)cos 5 - - (A . - G .)sin 5 . (9. 10) 


pxj xj pxj xj j pyj yj 


A^ . - G^. = (A^ . - G^.)sin 5 . cos b. + (A^ . - G^.)cos 5 . cos b. 
pyj yj pxj XJ J Pj pyj yj j ^j 


r r 

(A . - G Jsin B- 
pzj zj Pj 


(9. 11) 


A*^ . - G^. = (A^ . - G^.)sin 5 . sin b- + (A* . - G'‘.)cos 5 . sin b. 

pzj ZJ pxj XJ J Pj pyj yj j Pj 

(9. 12 ) 


+ (A^ . - G^.)cos B. 
Pzj ZJ J 


Substituting Eq. (9.7) into (9.9), writing the result as two scalar 
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equations (note that ^ . = P j Jj and solving for P . and R we 

cj ^ pj pj 

obtain 


P . = - Q . R . +7^ M . - p m(AP . - G^.)] 
PJ PJ PJ P»a c pzj zj 


(9. 13) 


R.=P.Q.+7^[M .-p m(AP.- G^.] 

PJ Pj PJ Ijj PZJ c pxj xj 


(9. 14) 


where M . and M .are the x . and z . components of external 
pxj pzj PJ PJ 

— th 

moment ^pj j blade, not including gravity. 

There is no equation for Q . because we have neglected the blade 

PJ 

inertia about the y . axis. Q the angular velocity of the blade -span 
PJ PJ 

axes along y depends explicitly on the blade pitch control input and is 
PJ 

not a state variable. It may very well be negligible. 

Equations (9- 13) and (9. 14) are integrated to obtain P . and R . , 

PJ PJ 

th 

respectively, for the j blade. They must be repeated N times for the 
N blades (j =1, 2, . . . , N). 


A. 2 . 1 0 A ERODY NAM IC FOR CES ON BL ADE S EGMEN TS 

In Section A. 2. 8 we developed the equations for determining the 

tangential (U^jj^), radial and perpendicular (Upjj^) components 

of velocity of the k segment of the j blade. These components 

are taken, respectively, along the blade span axes x . , y . and z ., as 
^ ^ e PJ "^PJ PJ 

th 

shown in Figure 10. 1. The total velocity of the k blade element with 


respect to the inertial reference frame is which is given by 

V = (U^ • + • )^^^ = |v 1 

^jk ^^TJk ^RJk ^Pjk^ ' ""jk* 


( 10 . 1 ) 


The yawed angle of attack between the segment velocity vector and 
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its projection on the blade -span -axis plane is denoted by 

in Figure 10, 1. Clearly 


sin <|> 


^_Pik _ ^Plk 

■ Iv.J - v.^ 


( 10 . 2 ) 


and 


(U^ + 

cos 4> = ^ 

jk 




( 10 . 2 ) 


Also, from Figure 10. 1, y , the angle -of -attack yaw angle in the x 

PJ 

y . plane, is given by 
PJ 


sin Y. 


rjk 


" (U^ , ,)i/2 

' Tjk Rjk' 


(10.3) 


or 


U 


cos Y.i 
jk 


Tjk 


( u ^ + tj 2 )^/2 

' Tjk ^Rjk' 


(10.4) 


th 

Now in general at the k segment, the blade itself has a pitch angle 
with respect to the blade -span axes, measured in the ^pj as 

shown in Figure 10. 1. From the Figure it is clear that 


tan 0 


(AB) _ (AB) 


Ajk (BC) (B*C)cos Yji^ 


(10.5) 


We now resolve 0 into the plane containing the yawed angle of attack 

./tI k 

9 . . This resolved angle 9 can be computed by noting that the dis- 

jk jk 

tance (A'B') = (AB) and 
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( 10 . 6 ) 


tan 6 - 


{A'B^ ) (AB) 


'Yjk (B*C) (B*C) 

But from Eq. (10.5), (AB) = (B*C) tan 0*., cos y-u* Substituting this 

Aj K j K 

into Eq. (10. 6), we have 


tan &Y = tan ^jk 


(10.7) 


The total yawed angle of attack, a , is simply 

jk 

O' = (b "b 0 

Y ^Y. Y 

jk jk jk 


( 10 . 8 ) 


or 


tan cj) + tan 0 

ik ik 

tan 

jk 1 - tan <p tan 0 

jk ^jk 


(10.9) 


where 


U 


tan (|) 


Pjk 


^jk 




( 10 . 10 ) 


or from Eq. (10. 4) 


U 


tan cj)^ 


jk 


Pjk 

iC7 ^jk 

T]k ■’ 


( 10 . 11 ) 


Substituting Eqs. (10.7) and (10.11) into (10.9) and solving for a^. , 

jk 

we obtain 


'Y. 


jk 


-1 ^Ajk^^Pik^l^°"Vlkl 

2 

U„., -U^., tan 6. cos v., 

Tjk Pjk Ajk 'jk 


( 10 . 12 ) 


where the absolute value signs are added around cos 1^ accordance 


60 



with the NASA Ames model. Also, the magnitude of the projection of 


Q>„ into the x z . plane is given by 

n PJ 

"TRANS., " '"y„ <=°“ 'jk' 

jk jk 

This X y . projection of angle of attack, namely 

PJ ^pj ^ ^ ^ TRANS., 

Jk 

as the angle of attack variable in the calculation of blade - section lift 


(10. 13) 
, is used 


and drag coefficients for small angles of attack. For larger angles a 
variation of Eq. (10. 13) is used, so that in general we write^ 


’'tRANS., Yji,.) 

Jk Jk ■' 


(10. 14) 


The lift and drag coefficients also depend on the Mach number, 
th th 

M., , for the k segment of the j blade, which is given by 


M.. 


(U^ + 

^ Rjk Pik^ 


(10. 15) 


where a is the speed of sound. Note that the Mach number is based on 
the velocity component in the ^pj» ^p. and not the total velocity. 

The segment lift and drag coefficients are given, respectively, by^ 


^fjk ^2^"^TRANS^j^. ^jk^ 


( 10 . 16 ) 


and 


"^djk = “Yjk ^3 ^“tRANS.j^. ^jk^ 


(10. 17) 


with a correction factor in tip loss in the last blade segment. 

It is assumed that the aerodynamic drag force acts in the 
direction opposite to the segment velocity V., , as shown in Figure 
10. 1, and that the lift force is perpendicular to and lies in the 


^ Mackie and Alderete, 0 £. cit. 
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( 10 . 18) 


plane of the yawed angle of attack a^. , In general 

jk 

= I p 'jk Sjk 


and 


°jk “ 2 ^ \ ^jk ^djk 


( 10 . 19 ) 


th 


where p is the air density and the area of the k blade segment. We 

can resolve L., and D., into force components F^., , F_,, and , along 
jk jk ^ Tjk Rjk Pjk ® 

the respective blade -span-axes x y . and z .. From Figure 10. 1, 

^ ^ PJ PJ PJ 


we see that 


F-k = ^^ik ‘*’y ■ °ik ^ ’‘'ik (10.20) 

^ ^ jk ■' jk ^ 


^Ejk = <^jk ■t'Y^ - °jk +Y > 

^Pjk = - <^jk <=°“ % + “jk ♦y > 


( 10 . 21 ) 
( 10 . 22 ) 


Finally, we need to write the formula for the geometric pitch 
angle, of the k^^ blade segment. This is given by^ 

®Ajk "" ®CUFF ■ “^S ■ ®is 


+ 6 , [y_, - e ’] - B- tan 5 + K or. + K 5? 

1 ^2k Pj 3 j j 


( 10 . 25 ) 


where the swash plate rotation, is a control input and p^, and 

5 . are problem state variables. The remainder of the parameters in 
Eq. (10.23) are constants, where only the y^j^ is dependent on the seg- 
ment index k. 

Equations (10. 1) through (10. 4)* (10* 12) and (10. 14) through (10.23) 
are used to compute the blade segment forces. They must be repeated 


Mackie and Alderete, op . cit. 
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Ns times for N blades with s segments per blade (k=l, 2, . . . j s). 

From Eqs. (10. 2) through (10. 4) we can rewrite Eqs. (10. 18) 
through (10.22) in the following form: 


^Tjk = if “Pik Vik- =difc“Tjk^ 


( 10 , 2h) 


*'Rjk = I \ ^jkK^ijk "Pjk ^jk - '"dik “Pik’ <“>• 

^PJk = b ^ VfJk'^Tik ♦'4jk'‘^"* VFjk’ 


These equations may be preferable to Eqs. (10. 18) through (10.22) 
for computer mechanization. 

Note that we have not included downwash corrections in the seg- 
ment velocity components ^Rjk’ ^Pjk respect to the 

air flow because of uncertainty on just how to incorporate the correc- 
tions. It appears that downwash effects, which certainly must be in- 
cluded, will not add a major computational load. 


A. 2.11 BLADE MOME N T EQUATIONS 

In the blade equations of motion in Section A, 2. 9 given by Eqs. 

(9. 13) and (9. 14 ) we see that the external moment components (not 

including gravity), M . and M . along the x . and z . blade span axes, 

^ ^ PXJ pzj ^ PJ pj 

must be calculated. In general the external moment vector M . consists 

PJ 

of aerodynamic and hinge moments, as expressed in the following equa- 
tion: 

M . = M + isl . (11. 1) 

pj aero ninge 
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Jir- 


where all moments are taken about the hinge point P. in Figure 1,1. 

th 

The aerodynamic moment due to the k blade segment is simply the 

cross product of the position vector, • from P. to segment c.p. 

2 k pj j 

th 

times the aerodynamic force at the k segment. The force components 
are given in Eqs. (10.20) through (10.22) or (10.24) through (10.26). 
The total moment is the sum of the segment moments. Thus 


M 


aero 


^ 2 k ^ k ^pj ^Rj k 


or 


M 


aero 


s 

^Zk^Pjk ^2k ^Tjk 




PJ 


( 11 . 2 ) 


The hinge moments are due to the spring-damper constraints for the 

lagging and flapping degrees of freedom of the blade. Thus M , the 

IjUj 

th 1 

hinge lagging moment for the j blade, is given by 


“LDi = + 


( 11 . 3 ) 


Reference to Figure 1, 1 shows that for positive 5 . the opposing 

moment due to spring restraint will act in the direction of and 

hence has a component M_ _ . cos B. along z .. 

^ LDj ^ PJ 


MpDj» hinge flapping moment, is given by 


(11. It) 


Reference to Figure 1. 1 shows that for positive the opposing moment 

lies along the x . axis. Thus 
PJ 

M,. = [K-0. + K-|S.]i . + L- [f ( 5 .) + sgn 5 .F'. ]cos 0.k . (11.5) 

hinge P'j P'j-’ PJ L' °j' ® J S ^J PJ 


Mackie and Alderete, _op. cit. 
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Combining Eqs. (11.2) and (11. 5), we obtain the following equations for 


the respective components of M along x . and z . : 
^ Pj PJ PJ 


^Zk^Pjk " ^FDj 


( 11 . 6 ) 


M . = V' 

pzj - ^ 


(11.7) 


where M, 


^ y^i E_., + M_ cos fi. 

^2k Tjk LDj 

and M . are given in Eqs. (11. 3) and (11. 4). Equations 


^FDj LDj 

(11.6) and (11.7) are repeated N times for N blades (j =1 , 2, . . . , N). 


A. 2.12 EQUATIONS FOR COMPUTING AND 

In Section A. 2. 2 we developed the equations for computing the 

components of angular velocity of the rotating shaft axes, x^^, 

y ,, z and hence the axes x' y* z* . in Figure 1. 1. Because ^2 

rj' rj ^ rj rj ^ ri 

th 

is the same for all N blades, the subscript j denoting the j blade is not 

included in Equations (9. 13) and (9- 114 -) in Section A. 2. 9, when 



integrated, give components of ^2 . , the angular velocity of the j 

^ ^ . • th 

blade -span axes. The angular rates 5^ and jS^ for the j -blade lagging 
and flapping angles are related to the difference between . and ^ 

1 . e # , 






- ^2 . 
pi. ri 
J 

.th 


( 12 . 1 ) 


where ^ is the cingular velocity of the j blade -span axes, x 

pr. ° PJ PJ 

^ th 

z with respect to the j rotating shaft axes, x y ., z or 
PJ ^ JO rj rj rj 

equivalentlvj x* y* .> ,. From Figure 1. 1 we see that 

^ ^ rj ^ rj rj 


65 


( 12 . 2 ) 


S2 = - 6 .k' . - S.i . 

J rj pj 

where k' . is the unit vector along the z' . axis and i . is the unit vector 


along the x . axis. 
PJ 


rj 


PJ 


Also, from Figure 1.1 it is apparent that 
ft* . = cos B. k . - sin B. j . 

PJ "^j "pj 

Substituting Eq. (12.3) into (12.2), we obtain 




jS. i . + 5. sin /3. j . - 5- cos S. k . 
J PJ J J “^PJ J ^J PJ 


Also, from Eq. (12. 1) we can write 


= (P . - pP.)i . + (Q . - qP.)J . + (R . - RP.)k 

Pi-j PJ rj PJ PI rj PJ PJ rj pj 

■ th 


(12.3) 


(12. U) 


(12. 5) 


where P ., Q ., R . are the j blade -span axis (P-frame) components 
PJ PJ PJ ^ 

of and P^., Q^., R^. are the j blade -span axis components of Q .. 

pi rj rj rj ^ ^ ^ ri 

Equating coefficients of i^^ and k^^ in Eqs. (12. 1^.) and (12. 5) and solving 

for 5 . and S., we obtain 
J J 


5 - ~ n • 

J cos jS. rj 


(rP. - R .) 

PJ 


( 12 . 6 ) 


and 


i3. - P^. 
J rj 


P . 
PJ 


(12.7) 


Equations (12. 6 ) and (12. 7) are integrated to obtain jS. and 5 ^., 


respectively. 


A. 2. 13 COMPUTATION OF HINGE SHEAR FORCE 

Each blade has acting on it the sum of the gravity force, m G, the 

aerodynamic force, F , and the hinge shear force at point P., - F . 
^ aero ^ J Pj 

The sum of these forces must equal the blade mass m times the 
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acceleration of the blade c. g. , A . . Thus, 

3 


mA . = mg + F - 

Cl. aero P. 

3 3 


(13. 1) 


where we have written the hinge force acting on the blades as -Fp 

3 

since the opposite force, Fp , will be the blade force acting on the hinge. 

- j 

Solving for Fp in Eq. (13. 1) we have 


Ft 3 = - m(A . - G) + F 

P. ci. aero 

3 3 


(13.2) 


In Section A. 2. 7 we have developed the equations for A . - G, the ac- 

celeration less gravity of the hinge point P.. From this we must now 
determine the acceleration less gravity of the blade c.g. By direct 
analogy with Eq. (5. 4) in Section A. 2. 5 

A. - G = A. - G+ ^ 2 . xp . + U . xp . (13.3) 

Cl. pi. pi. ^ cj pi. pi. cj 

3 ^ } ^3 ^3^3 

where p . = p j ., the position of the blade c.g. with respect to the 
CJ '^c'^pj ^ ^ ^ 

hinge point P.. Noting that Q . = P^. i . + Q^. j . + R^. k . and 

3 pij PJ PJ P3 PJ PJ P3 

U . = P^. i . + Q^. i . + R^. k we obtain from Eq. (13. 3), the 

Pij P3 P3 P3 P3 PJ PJ 

following_three scalar equations: 


cxj 

. 

cyj 

. 

czj 



= aP . . 

-gP. 

+ p 

io 

1 

in- 

(13. 

h) 

XJ 

pxj 

XJ 


c PJ PJ PJ 

-GP. : 

= aP . - 

-gP. 

- P 

2 2 
(R^. + P^.) 

(13. 

5) 

yj 

pyj 

yj 


c P3 P3 

- G^. : 

= aP . ■ 

-gP 

+ p 

(Q ,R . + P .) 

(13, 

6) 

ZJ 

pzj 

Z] 


c PJ PJ PJ 


'P -.P.. aP 

xj PV 
(9. 10) through (9. 12). 


where . - G^., A^ . - G^., and A*^ . - G^. are already given in £)qs. 
pxj xj pyj yj pzj zj 
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F is, of course, just the sum of the aerodynamic forces on 
aero 

each blade segment, as given in Eqs. (10.20) through (10.22) as com- 

th 

ponents k' k ^pj’ ^pj' ^pj ^ segment. 

Thus the hinge shear force components F F . along x , y , 

ir'Xj Jr'yj Jr'Zj p p 

z become, from Eq. (li. 2) 

P 



m(A 


cxj 


-gP.) 


s 

k=l 


Tjk 


(13.7) 


Pyj 


- m(A^ . - G^.) + Z) 


cyj yj 


k=l 


Rjk 


(13. 8 ) 


. = - m(A** . - qP.) + Z) F 

Pzj cz] zj P]k 


(13.9) 


These force components form the output of block 13 in Figure 2. 1. 


A. 2. U TOTAL ROTOR FORCE 

The hinge shear force components, F^ F^ and F^ for the 

Pxj Pyj Pzj 

th 

j blade, as given in Eqs. (13.7) through (13. 9) in the previous section, 
must be summed for all N blades to compute the total rotor force acting 
on the aircraft. But first the individual blade forces must be resolved 
into rotating shaft axes, i. e. , 

F^ = [T ]fP = [T (U. 1) 

P. pr. P. rp. P. 

3 3 3 3 3 

where [T ] is given in Eq. (2. 18). Thus 

F^ . = F^ .cos 6 . + F^ . sin 5 . cos + F^ . sin 5 . sin (Ik. 2 ) 

Pxj Pxj °j Pyj j "^J Pzj j '^j 

F!*^ . = - F^ . sin 5 . + F^ . cos g. cos B. + F^ . cos 5 . sin J3. (Ik. 3) 

yj xj J Pyj J ^j Pzj j ^j 
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'} ■ 

T 


. sin ^. + fP . cos )3. 

Pzj Pyj Pzj 


(U. It) 


Next the components must be resolved into shaft axes, which are com- 
mon for all blades, i, e. , 


= [T ]F^ = [T ]*^F^ 

P. rs.-* P. sr.-* P. 


(14.5) 


J 


J 3 


3 


3 


where [T ] is given in Eq. (2. 10). Thus 


,s 

= F^ . 

sin 

- F^ . cos 


(lit. 

6) 

Pxj 

PX] 

3 

Py3 

3 

s 

= F^ . 

cos 

+ F^ . sin 


(14. 

7) 

Pyj 

Pxj 

3 

Py3 

3 


,s 

Pzj 





(14. 

8) 


Summation over the N blades produces the total rotor force components 
along shaft axes. Thus 

N 

(U.9) 


f" = S f" 

Rx . , Pxj 
J=1 ^ 


N 

F® = S F® 
Ry j=i Pyj 


N 

F® r S F^ 
"'Rz "^Pzj 


(lit. 10) 


(lit. 11) 


To obtain the body axis force components, we use the transformation 

(lit. 12) 


F^ = - f^hs^X 


R 


where [T . ] = [T , ] because hinge and body axes have the same direc- 
sb s h 

tion. Using Eq. (2.3) for obtain 
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(U. 13) 


•—b . , _s . . ... _s 

F„ = F„ cos 1 - + F„ sin i , sin i„ + F.„ cos i, sin i_ 

Rx Rx 0 Ry 4> 0 Rz 0 


b s . „s 

F_ = F„ cos 1 - F sin i 

Ry Ry Rz (}) 


(14. U) 


^0 ^ 4y S ^0 ^Rz ^ci, ^0 

Equations (14.2) through (11^.. k) and (I 4 . 6) through (Ilf. 8) must be 
repeated N times for N blades. Equations (I 4 . 9) through (I 4 . 11) and 
(li|.. 13) through (14* 15) complete the rotor force computation. 


A. 2. 15 TOTAL ROTOR MOMENT 


th 


The moment transmitted through the j hinge for each blade con- 


sists of hinge spring damper components - M^^^^ along and M^^^. 

along X , in Figure 1.1, where M.^ and M.^.^. are given in Equations 
pj LDj FDj 

(11.3) and (11. 4 ). These moments are transmitted through to the hub 
along with the moment resulting from the hinge shear force calcu- 

lated in Section 13. Thus the total moment M, . acting at the hub for 
the j blade is given by 


^hj 


M_ k . + M^^y^.i . + e. X F_ 

LP] rj FDj PJ J P. 


(15. 1) 


where e^ = ej and is the vector position of hinge point with respect 

to the hub. We continue to neglect any monents along y' . and hence y .. 

The component of M_._ . i . along x* , and hence x . is simply M cos 5 . , 

FDj PJ ^ rj rj ^ ^ FDj ^3 

Equating scalar components along j and in Equation (15.1), 

we obtain 


r r 

M . = M cos 5 . + e F_ . 

hxj FDj Pzj 


T^. - e f!!, . 
hzj LDj Pxj 


(15.2) 

(15.3) 


70 



We resolve these components to shaft axes with the formula 


“h. = 

(15. 1^) 

j 

J J 


where [T ] is given in Eq. (2 
s r 

. 10). Thus, since we 

X 

assume M, . = 0, 
liyj 

<xj 

r 

= M, . sin 

nxj J 

(15.5) 

<yj 

r 

= M.-. cos 

hxj J 

(15. 6) 

<zj 

II 

(15.7) 


Since the shaft axes are common for all blades, the moment components 
in Eqs. (15. 5) through (15. 7) due to an individual blade can be summed 
over the N blades. Thus we have for the total rotor moment components 
at the hub the following equations: 


N 

s s 

J=1 


(15. 8) 


N 

s V' s 

1=1 


(15. 9) 


N 


M" - 


hz - 

1=1 


(15. 10) 


These are resolved into hub and hence body axes using the formula 


where is given in Eq. (2.3). Thus 

^hx ~ "^^x ^ ^6 ^ ^<j> ^6 


(15. 12) 
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(15. 13) 


M, = M, cos 1 . - M, sin i. 
iiy ny <p hz 9 


^z = * ^0 ^y ^<f) ^0 "■ ^z "e 

These hub moments are transmitted to the aircraft c.g. Also transmit- 
ted to the aircraft c.g. is the moment r x F , where F is the rotor 

Jrl XV K. 

force acting at the hub and computed in Section A. 2. 14 j and r is the 

ri 

moment arm from aircraft Cog. to hub. Thus the total rotor moment 

M acting about the aircraft c.g. is given by 
xv 


M 




(15. 15) 


where EqSo (15. 12) through (15. 14 ) give the components of hub moment 

h, , and F_ 


= F,_ i , 


along body axes, i ^ + y^^ j ^ ^ 

b — b — 

+ j b ^RzSd* ^^aking these substitutions into Eq. (15. 15), we 

obtain 


M 


M 


M 


Rx 

b 

Ry 

b 

Rz 


a , ^b „b 

^x ^ ^Rz - ^Ry 

(15. 16) 

Ky ^ - "h 4z 

(15. 17) 

^z ^Ry “ ^ Rx 

(15. 18) 


These are the rotor moment inputs to the rotational equations of motion 
in Section A. 2 . 1. 
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APPENDIX B 


SCALING CONSIDERATIONS IN SOLVING THE HELICOPTER 

FLIGHT EQUATIONS 

B.l, BLOCK 1 

In solving the helicopter flight equations using either a hybrid 
computer or a fixed-point digital computer, the scaling relationships 
between problem variables and computing -machine variables must be 
considered. Fortunately the maximum range of each problem variable 
in the equations of motion is either known or easily estimated. Under 
these conditions one can define new dimensionless problem variables 
given by the original variable divided by the maximum magnitude of the 
variable. Then by letting one computing machine unit equal one dimen- 
sionless problem unit, one can easily mechanize the equations. In the 
following sections we will discuss the application of this approach to each 
of the 15 computational blocks in Figure 2.1 in Appendix A containing 
the equations given in Sections A. 2. 2 through A. 2.15. 

First consider the 3 translational equations given in Section A. 2.1, 
Eqs. (1. 12i (1. 13), and (1.14). First we divide the equations by the gravity 
acceleration g^, so that every term in the equations is a dimensionless 


in units 

of go- 

Thus 

Zk 

W Q 
b'^b . 



+ 

So 

So 


- sin 0, + 


bx 


b Mg 


( 1 . 1 ) 


U,R, W, P, F, 

b b b b ... by 

+ + cos 0- sm \ 

go go b Mgo 


( 1 . 2 ) 
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b b b , b b , „ bz m 

= - + + cos 0, cos 4>, + ^ (1. 3) 

g g g b b Mg^ 

®o ®o ®o o 

Rather than use a separate maximum range for each variable, it 
would appear simpler to use a single maximum range for each group of 
three variables representing scalar components of a single vector. For 
example, if maximum expected magnitude of the vehicle 

velocity vector, then we define dimensionless body-axis velocity com- 
ponents (using lower case letters) by the following formulas; 


U, 





(1.4) 


max max max 

Similarly, for maximum magnitude of the body- 

axis angular velocity vector, we obtain the following formulas for dimen- 
sionless angular velocity components: 



max max max 


Although this scheme is not quite as efficient as introducing a separate 
maximum value for each individual component, it is simple and should 
easily provide sufficient accuracy when a 16 bit digital word is employed 
for each dimensionless variable. 

Let us also define dimensionless external body-axis force com- 
ponents as follows: 


^x 



max 




Mg., 




t = 


bz 


max 


bz Mg 


'S%>, 


»max 
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Using Eqs. (1. 4) - (1. 6), Eqs. (1. 1) - (1. 3) become 

g 


''b'b-‘'g ="®b * 


t -Vb+ '"bPb^ ■'g ■=ose^sm+|^+ 

'"b'“ma, =t ■''b‘’b'- Vb-" *=g ““®b 

where the constants and are given by 

k = 


(1-7) 
( 1 . 8 ) 
(1.9) 


g 


g F 

2 k = k 

max max ^ o max 


( 1 . 10 ) 


and will both be less than unity for any practical case. 

All the variables on the right side of Eqs, (1. 7) thru (1.10) are 
dimensionless and easily interpreted. To retrive the 

scaled variables (e. g. , for display purposes) one need only multiply by 

the maximum value used to make each variable dimensionless. The 

common factor Eqs. (1. 7) - (1. 9) is simply the integrator gain in 

an analog mechanization. Jh digital integration Is multiplied by the 

numerical interval of integration T, where clearly the product ^ ^ 

in any practical case. It also seems clear that the sum of the bracketed 

terms on the right side of Eqs. (1-7) - (1. 9) is quite unlikely to exceed 

unity. This is because terms such as Eqs. (1.9) have a maximum 

possible value of unity, which represents an actual acceleration equal to 

V ^ . Typically, we might have V = 300 ft/sec, ^ =1 rad/ 

max max ^ ^ max max 

2 

sec in a helicopter simulation. Then V = 300 ft/sec or almost 

^ max max 

10 g* s . Even if both \a^ and q^ reached maximum values simultaneously, 
i. e. , corresponding to an actual acceleration magnitude of 10 g’ s, the 
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the remaining terms within the brackets in Eqs. (1* 7) - (1. 9) woiild have 
much smaller equivalent accelerations in any allowable flight maneuver. 


Next, we consider the rotational equations of motion given by 
(1. 2Z), (1,23), and (1.24) in Section A. 2. 1. We define dimensionless 
external moment components by 



( 1 . 11 ) 


Using Eqs. (1, 5) and (1.11) above the rotational equations (1.22), (1.23), 
and (1. 24) in Section A. 2. 1 become 


^2 


( 


max 


I -I 

zz yy V 

XX 




max 


I -I 

, XX zz , 

= - ' “i ' 'bPb ''m“b, 

yy 


( 1 . 12 ) 


(1.13) 


where 




max 


, yy XX , 


zz 




max ^ 


max 


(1. 14) 


(1.15) 


All terms in Eqs. (1.12, (1.13), (1.14) are dimensionless and will 
remain less than unity, with the ^m^^y* ^m^^z 

predominating (the effect of external moments), p^, q^, and r^ are 
obtained by integration of the respective equations using an integration 


gain of ^2 


max 
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In terms of dimensionless angular rates the Euler angle equations 


(1. 2 5), (1. 26 ), and (1. 27 ) in Section A. 2. 1 become 






max 


= (r, cos(|)_ + q_ sin<j), )/ cos0 

D D b b b 


(1.16) 




max 


= cos<t>^-r^ sini^^ 


(1. 17) 


= Pb + n ®“®b 


(1.18) 


max 


max 


Again, integration with gain yields 0^, (|>^. Actually, the gain 

will be ^ /tv so that the resulting ajigles range over + 1, correspond- 
ing to IT radians. 

Finally, the equations for earth-axis velocities, (1,30), (1,31), and 
(1. 32) in Section A. 2. 1, become 
S 


V 


max 


u^cos 0^cosi|j^ + v.^( -cos(j)^sin4j^+ sinc|)^sin0^cosilj ^) 


+ w^(sinc|i^ sini|j^ + coscfi^sin 0^cosi|j ^) 


(1- 19) 


JL 
T 

max 


= Uj^cos 0^smi|j^ + Vj^(cosc|jj^cosi1j^+ sin<ji^sin0^sinilj ^) 

+ w^( -sinc|)^cosiIj^+ cos(|)^sin0^sinilj^) 


( 1 . 20 ) 


— = = -u, sm0, + V, sind), cos 0- +w, cos 0, cos4>, 

V V ^ b b ^b b b b ^b 

max max 


( 1 - 21 ) 


Here the integrator gains will depend on the desired range in distance 

north (S ), distance east (S ) and altitude (h = -S ), 

X y z 
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B, 2 Block 2 


Here we compute the angular velocity components of the rotating 
shaft axes. As before we use lower case letters to denote dimension- 
less components. Thus 



max max max 


Then Eqs. (2. 7), (2. 8), and (2. 9) in Section A. 2. 2 become 


g 

p = p, cos i - r, sin i (2.2) 

r 0 b e 

q = p, sin i , sin i^ + q, cos i . + r^ sin i , cos i^ (2. 3) 

^r ^b (|) 0 ^b ct> b f|) 0 

n 

r^ - ■ ^ = p^ cos i^ sin i^ -q sin i^ + r_ cos i , cos i _ (2.4) 

r^^ ^bd) 0^b<bbcb 0 

max 

s s 

Note that r is computed rather than r . Thus we compute the 

r s max r 

dimensionless difference between the z axis components of rotating 

s 

shaft axis angular velocity and angular velocity of rotating shaft axes 

Q 

with respect to shaft (non -rotating) axes. In this way r does 

^ ® ' r s max 

g 

not exceed unity, whereas r would, since {^ ) > > ^ 

^ r s max max 

Eqs. (2. 12), (2. 13), and (2. 14) in Section A. 2. 2 are simply re- 
written in dimensionless form, i. e., with lower-case letters. Thus for 
the jth blade 


17 S S 

p = p siniLi . + q cosdj , 


r s s . 

q . = -p cosijj . + q sim|j . 
Tj r J r j 

r s _ s s 
^ max ^ max 


(2. 5) 
( 2 . 6 ) 
(2.7) 
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Finally, from Eqs. (2. 20) and (2. 22) in Section A, 2. 2 we have 
the following expressions for dimensionles s .blade -span axis components 
of rotating shaft-axis angular velocity; 


r r 

p‘ . = p . cos 6 .-q . sm 6 , 
rj ^rj j ^rj j 


( 2 . 8 ) 


12 

.P L, 


max 


r r A 

r‘ . -t; cosp . = p . sin 6 .sinp . + q . cos 6 .sinp .+ (r .- 7 ^ )cos p 

ri n J rj j "^3 3 ^3 rj 


n 

r s 


max 


(2.9) 

where we have not included the formula for q^ since it is not needed. 

B. 3. Block 3 

In this block the angular acceleration components of the rotating 
shaft axes are computed. In terms of the dimensionless angular accelera- 
tion components given in Eqs. (1,12), (1.13), and (1.14) in Section B.l 
Eqs. ( 3 , 6 ), (3,7), and ( 3. 8 ) in Section A. 2. 3 become 


Ps Pb . ^b 

* cos 1 


12 12 
max max 


B n ^0 

max 


(3.1) 


^s Pb . . . . % 

smi^smig + — 


12 12 
max max 


cos i . + 


max 


4 


max 


sm 1 , cos 1 ^ 
(j) 0 


(3.2) 


. s 

r 

s 


Pb % ^*b 

n S'"”® 'e 


max max 


max 


max 


(3.3) 


Eqs. (3.17) and ( 3. 19) in Section A, 2 . 3 for the jth blade become 




12 


max 


. s 

^ s 

^ sindr . + cosd . 

12 1 12 1 
max max 


• s 

Ps 


(3.4) 
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-n /S2 

rj s max 


(3. 5) 


. s 

r 

s 


Q 




max 


max 


Note that r , -^2 /^2 is computed rather than r . . This is because 
rj s max rj 

in general (^ ) » ^ , i. e. , the angular acceleration of the 

^ s max max 

rotating shaft axes with respect to the shaft axes (non -rotating) can be 

much larger that the overall aircraft angular accelerations. Thus 

T* * r 

r -^2 /^2 will remain reasonable in magnitude, whereas r . by itself 
rj s max rj 

will not. 

B, 4, Block 4 

In Block 4 of Figure 2. 1, Section A, 2. 2, the hub-axis velocity 
is computed by adding the velocity of the hub relative to the body axis to 
the body axis velocity. The relative velocity components are normalized 
by dividing the actual components by maximum value of the 


hub -to -body relative velocity. Thus 


b ^hb 


hb 


W, 


'\b " V * ■^hb ~ V ’ ^hb " V 

r max r max r max 


hb 


(4,1) 


If d^^^ is the largest of hub relative position components 
then clearly 


V = ^2 d 

rmax max max 


(4.2) 


In general V will be small compared with V , the maximum air- 

® rmax max 

craft velocity. From Eqs. (4.1) and (4.2) along with Eqs. (1.4) and (1. 5) 
in Section B.l, Eqs* (4,4), (4. 5), and (4.6) in Section A. 2. 4 for the hub- 
to-body axis relative velocity components become 
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b 



\b 

d 

max 

\"d ’^b 

max 

b 



'"hb 

“d 

max 

’^b ■ d % 

max 

b 



^hb 

d ^b " d % 

max max 

(4.9). 

(4.10). 

and (4.11) in Section 


(4. 3) 


(4. 4) 


(4. 5) 


following for the dimensionless hub-axis velocities along the shaft axes: 

(4.6) 




(4.7) 


w 


^ b b b 

*^=(u, +k u, _ )cos i, sin -(v, +k v, - )sin + (w- +k w. _ )cos i . cos 
r D r nb cj) 0 b r hb c|) ' b r hb ^ 0 

(4,8) 


where 


A 

k = 
r V 


rmax 


max 


Q d 
max max 

V 

max 


(4.9) 


s s s 

In Eqs, (4. 6)-(4. 8) we have used V to make u , v , and w dimen- 
^ max r r r 

sionless, i. e. , 


u® 

u® =-^ 
r V 


max 


V 

s r 

. V = — — 
r V 




w 


max 


max 


(4, 10) 


Note that the constant k^ will generally be small compared with 


unity. 

B. 5 Block 5 

in this block the hxib-axis accelerations are computed from the 
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body-axis accelerations less gravity. The latter are simply the external 

forces acting on the aircraft divided by the aircraft mass. Thus the 

maximum magnitude of body-axis acceleration less gravity is simply 

(F, /M) . This is then used to define dimensionless components of 

b max 


body-axis acceleration less gravity, i. e. 




A, -G 
bz z 


A, -G 

by y _ 

\x“®x ■ (F /M) * ^y'®y"(F, /M) (F /M) 

b max b max b max 


We also use (F. /M) 

b max 


(5.1) 

to define dimensionless components of hub-axis 


acceleration less gravity. Thus 


A -G 
hz z 


A^ -G A^ -G 

hx X hy y _ 

^hx“®x "(F ,/M) »®’hy"®y"(F /M) ’^z"®z"(F /M) 

b max b max b max 


(5.2) 


From Eqs, (1. 5) in Section B. 1, Eqs. (5.1) and (5.2) above, and Eqs. 
(4.3), (4.4), and (4. 5) in Section B. 4, Eqs, (5.5), (5.6), and (5.7) in 
Section A. 2. 5 become 

. b 


a, -g = 


b 




^b 

. \ 

‘^b^hb 

■^d 


■ d 


) 


max 

max 

max 

max 

b 




Pb 

_ \ 

■Pb^hb 

■^d 


’ d 


) 


max 

max 

max 

max 

b , 


Pb 


% 

. \ 

■%“hb " d 

U 

■ d 


) 


max 

max 

max 

max 


(5. 3) 


(5. 4) 


(5. 5) 


where 


A ^ d 
u - max max 

a “ (F, /M) 


(5.6) 


b 'max 

For the dimensionless components of rotating shaft-axis accelera- 
tion less gravity we have from Eqs. (5.10), (5,11), and (5.12) in Section 
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A, 2. 5 


a® -g® = ( a, -g )cos i- - (a, -g ) sin i 
rx ®x ' hx ®x' 0 ' hz ^z' ( 


(5.7) 




a® -g®=(a, -g )cos i, sin i^-(a, -g )sin i +(a -g )cos i, cos i (5.9) 
rz z hx ^ e ' hy <j) ' hz ®z ^ 0 


B. 6. Block 6 

In block 6 the velocity of the blade -span axes, i. e. , the hinge 

point Pj for the jth blade, is computed for each of N blades (j = 1, 2, . . . , N). 

Again we use V (the maximum body -axis velocity) to normalize the 
® max 

velocity components, which in dimensionless form are given as lower- 
case letters. Also, we recall that ^ (the maximum body-axis angu- 

max 

lar velocity) is used to normalize the angular velocity components. Thus 
fromEqs. (6.2), (6.3), (6.4), (6.8), (6.9), and (6. 10) in Section A. 2. 6 
we obtain for the jth blade 


Q 

r s . . . s . , r s , max 

u ,=u smib .+v cosib .-(r .--rr ) t;: 

PJ r j r J ^ V 

J J max max 


e^^Q 


max 


max 


CO 


( 6 . 1 ) 


r s s 

V . = -u cosili . V sindj . 
PJ r r ^3 

r s . max r 

w . = w + p . 

PJ r V ^rj 

max 


( 6 . 2 ) 

(6.3) 


where the dimensionless angular velocity of the rotating shaft axes 


relative to the non-rotating shaft axes is given by 

^2 


CO = 


s ^ 


(6.4) 


max 
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Note that , as computed in Eq, (2. 7) of Section B. 2, 

rj s max 

is used in Eq. (6.1) rather than r^ . itself to preserve good scaling. For 
e, the distance from hub to hinge, equal to perhaps one or two feet and 


equal to perhaps 30 rad/sec, clearly e ^ /V will gen- 

s s max 

max max 

e rally be small compared with unity. Since » clearly 


n>ax 

e ^ /V is quite small compared with unity. Thus good scaling of 
max max 

Eq«s. (6.1), (6.2), and (6.3) is assured. 

B.7 Block 7 

Here we compute the dimensionless components of acceleration 
less gravity for the hinge point of each blade. First we write the di- 
mensionless components of rotating shaft axis acceleration less gravity. 
In Eqs. (5.7), (5.8), and (5.9) these are given as components along the 

shaft axes. Using the transformation matrix [ T ] in Eq, (2.10), Sec- 

s r 

tion A, 2. 2, these become the following components along the rotating 
shaft axes for the jth blade: 


r r ,s s,.,,,s s 

. -g . = (a- -g )smib.+(a -g ) cosil. 
rxj ^xj rx ^x j t-y ^y 

•(7.1) 

r r /S s. ,.s s.., 

• -g •= -(a- -g )cosi|j .+ a -g ) sm4j . 
ryj yj rx ^y' ry '^’y' 

(7.2) 

r r s s 

.-g = a -g 
rzj °zj rz 

(7. 3) 


Next Eqs. (7, 5), (7.6), and (7.7) in Section A, 2. 7 are used as the 
basis for computing the dimensionless components of blade -span-axis 
acceleration less gravity. We recall that normalizing 

factor for acceleration components and ^ for angular acceleration 
and velocity components. Also, Eq. (2. 7) in Section B. 2 shows that 


84 



r computed rather than r^ . directly. Similarly, Eqs, (3, 5) 


shows that (ir^ ~ Q /Q )/^ is computed rather than r^. 

rj s max max rj max 


We 


also note that co = /SI 
s s s 


in Eq, (6. 4), Section B,6, With these 


max 


considerations we obtain the following equations as the dimensionless 
equivalents to Eqs. (7, 5), (7, 6), and (7,7) in Section A. 2. 7: 


r r r r , , r r r 
a .-g . = a .-g .+ k [ q .p . - 
pxj ^xj rxj ®xj q rj^rj 


r' -n /£2 

r j s max 


o> 




max 


- ] -k k ^ 
■* q r 


(7.4) 


max 


.-g^.+k k^co^=a’^ .-g’^.-k [ (r^.-n /J2 )^+(p^)^] -2k k (r’^.-a /S2 Vw 

pyj yj q r s ryj ®yj q*- ' r j s max' •' q r' rj s max' f 


. r 

Prj 

a^ .-g^.=a^ .-g^.+k [ .(r^ )+ — ] +k k q^ .<o 

pzj ^Z3 zxj ®zj q^ ^rj rj s max q r rj s 


(7. 5) 
(7.6) 


where 


( S2 ) e 
max^ 

■q" (F^/M) 

^ ' b 'max 


Q 


k = 


* k = 
r 


max 


(7.7) 


max 


Typical parameter values might be 1 rad/sec, e = 2 feet, 

(F /M) = 100 ft/sec^ (~3g) , « = 30 rad/sec. Then k =0.02, 

' b max ' ® ’ s q 

2 2 

k = 30. Thus the term k k co in Eq. (7. 5) can range up to 18, so that 
r q r s 

subtraction of this centrifugal acceleration from the y component of 


hinge acceleration -less -gravity improves the scaling considerably. 

B. 8. Block 8 

In block 8 Eqs. (8. 5), (8. 6), and (8. 7) in Section A. 2. 8 are im- 
plemented to transform the blade -span axis velocity components from 
rotating shaft axes (R -frame) to blade-span axes (P frame). In dimension- 
less form the equations become 


85 



( 8 , 1 ) 


P c ^ • C 

u. . = u . COSO .-V . smo 


PJ PJ 


J PJ 


J 


= u^.sinS .cosB .+v*^.cos6 .cosB .-w^.sinB . 
pj PJ J J PJ J *^J PJ J 

w^.= u^.sinS ,sin3 ,+v^.cos 6 .sinB ,+w^.cosB . 
PJ PJ J J PJ J J PJ J 


( 8 , 2 ) 

(8,3) 


Eqs, (8, 8), (8, 9), and (8,10) in Section A, 2, 8 give the tangential, radial, 
and perpendicular components of the center of pressure of the kth segment 
of the jth blade, the tangential component, will in general be the 

largest, with a maximum value well in excess of the maximum 

aircraft velocity. Thus we use a normalizing factor for the 

blade-segment velocity components, where ~ 1000 ft/sec. Then 

the dimensionless components for the kth segment of the jth blade become 


^Tik 




jk 


u. 


u 


Pjk 


‘Tjk ’ "Rjk V, ’ “pjk V, 

bmax bmax ^ bmax 


(8.4) 


where these components are, respectively, along the x .,y .,z . blade- 

PJ PJ PJ 

span axes. 

Let ^ be the maximum value of the components P . , Q , , and 

pmax PJ PJ 

R .-^2 of blade -span axis (P-frame) angular velocity. We use ^ as 

PJ s pmax 

the normalizing factor for the dimensionless blade-span axis angular 
velocity-components. Thus for the jth blade 


-EJ_ 


Q . 
_Ei. 


^pj ^ ’ ^pj ^ 

pmax pmax 


, r .-k CO - 
PJ s s 


R .-^2 
, .P J ■ s 






smax 


pmax 


s ^ 


[S. 5) 


pmax 


where co = as defined in Eqs. (6.4), Section B.6. 

s s smax 

Noting that V is the normalizing factor for the velocity corn- 
max 


ponents in Eqs. (8.1), (8.2), and (8. 3), we have from Eqs. (8.4) and 
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(8^ 5) the following forncuilas equivalent to Eqs. (8. 8), (8. 9), and (8.10) 
in Section A. 2. 8: 


V ^ y 1 Q y 

max p pmax 2k , , , smax 2k 

"Tjk= T “pj - 


bmax 


bmax 


max p 


Pjk 


max 


bmax 


w^. + 

PJ 




smax 


bmax 


PJ 


( 8 . 6 ) 

(8.7) 


( 8 . 8 ) 


Here the coefficient V /Y. will be leas than unity whereas 

max bmax •' 

will be near unity for the blade tip segments. 

smax D max 

B. 9. Block 9 


In this block the blade equations of motion are represented. Eqs. 
(9.10), (9.11), and (9.12) in Section A. 2, 9 simply transform the components 
of blade -span axis acceleration less gravity from rotating shaft axes to 


blade -span axes. Rewriting the equations in dimensionless form using 
Eqs. (7. 4), (7.5), and (7. 6) in Section B. 7, we obtain 


P P /T r, - ,r r,,,2 2..-,,.22.^- 

a .-g . - (a .-g ,)cos6 .-(a .-g ,+k k co ) smo ,+k k co sin6 . 

pxj ^XJ pxj ^XJ J pyj yj q r s j q r s j 


(9.1) 


pp22 rr rr22 

.-g''^.+k k w cos6 .cosB .=(a ,-g .)sin6 .cosB .+(a ,-g .+k k w )cos6.cos6 

pyj yj q r s j pxj ®XJ J pyj ^yj q r s 3 


r r 

-(a .-g .)sinp . 

pzj ^zy ^2 

.-g^.= (a^ .-g^.)sin6 .sin(3 .+(a^ .-g^.+k k^ w^)cos 6 .sin(3 . 
pzj ^Z2 pxj ^XJ J J pyj ^yj q r s j ^j 

2 2 

+(a^ .-g^ .)cos P . -k k o) cos6.sinp. 
pzj zj ^j q r s J J 


(9.2) 


(9.3) 


where co = 12 , as given in Eq. (6.4), Section B.6. 

s s smax 
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.if 


M . and M . in Section A. 9 represent the x . and z . components 
pxj pzj ^ pj PJ ^ 

of external moment less gravity acting on the jth blade , respectively. We 
make them dimensionless by dividing by 


M 


m 


pxi 




M 


_E22 




(9.4) 


max 


pxj= (Mp/Ij^)n^ax ’ 

Then the dimensionless counterpart of the equations of motion for the jth 
blade given in (9-13) and (9-14) of Section A. 2. 9 become, with the help of 
Eqs. (8. 5) in Section B. 2, the following: 


J2L 


12 


smax 


q .(r .-k CO )-q .co +k m -kja^ 
c PJ PJ s s ^pj s m pxj f pzj zj' 


( 9 . 5 ) 


r .-k 
P.l s s 

12 ' k 

smax £ 

where 


p .q .+k m J-k-(a^ .-g^-)-7T 
PJ PJ na pzj r P3CJ ®xj' 


pmax 




max 


k = 
m 12 12 

pmax smax 


p m(F /M) 


S2 


max 




pmax smax 




(9.6) 


(9.7) 


Eqs. (9. 5) and (9. 6) are integrated with gain 12 to obtain p . and 

^ ^ ^ amax '^pj 

r^j-k^to^. They are repeated N times for N blades (j=l,2,. . . ,N). 


B.IO. Block 10 


In block 10 in Figure 2. 1 in Section A. 2. 0 the aerodynamic forces 
acting on the individual blade segments are computed using Eqs. (10.18) 
thru (10. 22) in Section A. 2, 10. First we define dimensionless tangential, 
radial, and perpendicular velocity components for the kth segment of the 
jth blade using the following formulas: 
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^Tjk ” J, 


Tik 


F. 




R.ik 


A -,r2 ’ Rjk 1 . 2 

^pA VT ^T_ 

2 max bmax 2 max bmax 




Pjk" 1 


F 

Pjk 


TP A V“ 

2 max bmax 


Assuming | C 


( 10 - 1 ) 

= 1, Eq, (10.1) wiU define dimensionless force 


i jk max 

components with a maximum value of unity. The normalizing factor 


can easily be adjusted to accommodate the actual maximum and 
range. From Eq, (10.1) above the dimensionless force equations can be 


written from Eqs, (10. 18) thru (10. 22) in Section A. 2. 10. Thus 


( 10 , 2 ) 


^ vf, [ C^,, sinct>v -C,., cos c(>„ ] sinv 


Rjk A jk^ ijk“ ^Y., djk ^Y., 

max jk jk 


jk 


(10.3) 


^jk’ 


max jk ■' jk 


(10.4) 


where 


,2 , 2,2 , 1/2 

'"jk ^'^Tjk '^jk'^ ^Pjk^ 


(10. 5) 


The formulas for cosy , siny , cost|) , sincj) remain as before ex- 

^ Jk 

cept for the substitution of dimensionless velocity components u » 

ijK Kjk 

u_., for , U . 

Pjk Tjk Rjk Pjk 

The formula for is reasonably scaled as is, with all angles 

in radians. The equation is rewritten using dimensionless velocities 

jk 

'^Tjk’ '^jk’ '^Pjk’ formula for Mach number at the kth segment. 

Thus 


89 



Ill lll■■ll I nil 


,, bmax ,2 .2 »l/2 „„ , 

“j ; <"Tjk + “Pjk> ' 

An alternative mechanization to Eqs. (10.2), (10.3), and (10.4) 

is based on Eqs. (10. 24) , (10. 25) , and (10. 26) in Section A. 2. 10. Thus 


^Tjk A ^ ^jk ^ jk *^jk ‘^°®'^jk"*^djk'^Tjk ^ 


^Rjk= A V jk'Vjk"“^jk-^djk'^jk ^ 


^Pjk A ^jkl^ jk^'^Tjk “Rjk^ ■'■^djk^I^k^ 
max 

Slight variations in the scaling suggested above may be needed to obtain 
near optimal scaling, but the basic approach would be identical. 

B.ll. Block 11 

Using Eqs. (9. 4) and (10. 1) we can rewrite Eqs. (11. 3), (11. 4), 
(11. 6), and (11. 7) in Section A. 2. 11 as the following two equations: 


pxj yk=iY2^^^ 


^ ^jk M 


K . n (3 . 

P 6 + P smax J 

[ M n 

pmax pmax s max 


(11. 1) 


s y. 


. = -k Z/ ( 


) frr,.-.. 


pzj 7k=l 72^^^ Tjk M 


pmax 


[ f(6^.)+sgn6^ F'^ ] 


- A 

2 max b max ^2 max 
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B.12. Block 12 


Next we consider the dimensionless form of Eqs, (12. 6) and (12. 7) 
* • 

in Section A. 2. 12 for 6 .and p.. From Eq. (8, 5) in Section B. 8 we note 

J J 

that ^ is the normalizing factor for P , and R , and from Eqs. 

p max PJ PJ s 

(2.1) and (2. 9) in Section B. 2 is the normalizing factor for and 

' max rj 

R^. cos6 .. Thus Eqs. (12. 6) and (12. 7) in Section A. 2.12 become 
r j s j 'i V 


1 ^ 

= —L—r-l (rP . 




^2 cosB . k r . 

pmax 3 r J 


max 


cosp .)-(r .-k <0 )-k co (1-cosp .)] (12.1) 

j pj s s s s j 


and 


p. k 

s p 

^ “ k ^rj " ^pj 

p max J, J 


( 12 . 2 ) 


Here ^ =^2/^2 (Eq. 6. 4 in Section B. 6), k = ^2 /^2 (Eq. 

s s s max r s max max 

7. 7 in Section B. 7) and k = ^2 /^2 (Eq. 8. 5 in Section B. 8). 

s s max p max 

We recall in Eq. (12.2) that the two z^^-axis dimensionless angular rates, 

r^.-(^2 /J2 ) cos P . for the rotating shaft axes (R-frame) and r ,-k co 

rj s max ^ j ^ pj s s 

for the blade span axes (P-frame), do not include the component due to 

rotor-shaft spin rate ^2 =co ^2 . I, e. , it is subtracted off of each 

3 8 s max 

component as computed respectively in Eqs. (2. 9) in Section B. 2 and 
(9- 6) in Section B. 9- This avoids computing the small difference of 


large quajitities to obtain 6y Equations (12.1) and (12.2) are integrated 

with gain ^2 to obtain P . and 6 .. 

pmax J 

B.13. Block 13 

To compute the dimensionless form of the acceleration compo- 
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nents of the c, g. of the jth blade, we use 2p ^ ^2 as a 

c s max p max 

normalizing factor. Recalling that (F /M) is the normalizing factor 

D max 

for acceleration components of the hinge point P. (see Eqs. 9.1|9.2,and 

J 

9. 3, Section B. 9), ^ is the normalizing factor for blade angular 

p max 

velocity components (Eq. 8, 5, Section B, 8), and normalizing 

factor for blade angular acceleration components (Eqs. 9. 5 and 9-6, 
Section B. 9), we can rewrite Eqs. (13. 4), (13. 5) and (13. 6) in Section 
A. 2. 13 as follows: 

P _„P ^ t /.p _„P r -L 


CO 


s s max pmax 


(13.1) 


.P 

cyj ''yj 


a‘ .-g^. + “■ (1 + ) CO = k [ a^ .-g^,+k k co cos 6 .cos p.] 


2 2 

t w 

pyj “yj q r s 
.2 


-T~ [ p^.+(r .-k CO )^] - CO (r .-k co ) + — - — ^ ( — )(l-cos6 .cos | 

2k_ pj pj s s s pj s s 2 p j 


2 .f I »2i 

3 ,+(r .-k CO ) 

PJ PJ s s 


k CO* 


aP .-gP.= k (aP .-gP.) +'7 [ 7^ q . (r .-k w ) + q .w +„ 

czj ^zj c pzj ®zj 2 k ^pj PJ s s ^Di s n 




PJ 


s max 


J 

(13.2) 

(13.3) 


where 


k =, 


(F, /M) 
b max 


c 2 p ^ ^2 

c s max p max 


(13.4) 


kg e 2 

In Eq. (13. 2) the term — (1 + — ) co represents the negative of 

2 p s 

^ c 

the centrifugal acceleration of the blade c, g. By substracting it from 


the total y^j axis acceleration component we improve the scaling by a 
very sizeable factor. 

For the normalizing factor for the components of hinge shear force 

let us use 2mp o2 ^2 . Thus 

c s max p max 
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rP ^ Pn f P „ PY J 

Pxj 2m p ^2 ' Pyj 2m p ^ ^2 

^ c s max p max c s max p max 


fP ^ ^£5J n 

Pzj 2m p ^2 12 ' 

c s max p max 

Using Eqs. (13.1) thru (13. 5) above, Eqs. (13. 7), (13. 8), and (13.9) in 
Section A. 13 become the following! 


fP . = .(aP ..gP.) + E k 

Pxj cxj ®xj' a Tjk 


• - kco^ = -(a^ ^ w^)+ E k f_ 

T^j s s ' cyj ®yj s s' ^=1 ^ 


f . = -(a'^ + Ij k^fo.. 

Tzj czj ®zj a'Pjk 


where 


1 A 

*:r p A V- 

2 ^ max bmax 

a 2mp 12 12 

^ c s max p max 


^ c 


Again note in Eq. (13. 7) that the hinge force component has 

' 2 

the centrifugal force berm k subtracted off. 

s s 

B.14. Block 14 

The dimensionless hinge force components for the jth blade, as 

given in Eqs. (13. 6), (13. 7), and (13. 8) in Section B. 13 must be resolved 

from blade -span axes (P-frame) to rotating shaft axis (R-frame) compo- 

nents using the matrix [ T ] , where [ T ] is given in Eq. (2.18), 

rPj Tp. 

Section A. 2. 2. Thus we obtain 


.= . cos6 .+(fl .-k w )sin6 .cosp .+fE .sin5 .sinp .+k w sinS.cosp. 

Pxj Pxj j ' Pyj s s' j "^j TPzj j j s s ,j '^j 


(14. 1) 
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j=-f^jSin5j+(f^yj-kgUg)cos5 xosp j+f^jCos6jSmp j 

(14. Z) 

.=-(f^ ,-k* )sinp .+f^ .cosp.-k* u^sinB. (14.3) 

Pzj 'Pyj s s' *^3 pzj ^3 s s *^3 

We note that k* as defined in Eq. (13. 9) of Section B. 13 is equal approxi- 
s 

mately to k /2 which from Eq. (8. 5) in Section B. 8 is equal to 
s 

^ /Z^ . We recall that ^ is the maximum rotor spin 

s max p max s max 

velocity and maximum angular velocity of the rotor-blade 

components after ^ has been subtracted. Thus ^ /2^ (i. e., 

s s max p max 

I 

k ) will be likely to exceed unity by an appreciable amount. On the other 
s 

hand k* only enters Eqs. (14.1) and (14. 3) when multiplied by sin 6 . or 
® J 

sinp so that the product should remain less than xmity. Of course the 

t 2 

k CO cos 6 ,cos3 . on the left side of Eq. (14. 2) will be much larger than 
s s j j ^ ^ 

unity. This is exactly why we compute directly the small difference 

* 2 

f .-k CO cos6 .cosp .. 
pyj s s j j 

The dimensionless force components given in Eqs. (14.1), (14.2), 
and (14. 3) are resolved from rotating shaft axes (R-frame) to shaft axes 


» 2 . 


(S -frame) using the transformation matrix [ T 1 , where [ T 1 is 

^ *■ sr. '■ sr. 

J J 

given in Eq. ( 2 . 10 ), Section A. 2 . 10 . Thus we have 

s t 2 r r > 2 * 2 

£Z .+k CO cosi]j .=f .sinijj .-(fp .-k co cos 6 .cosp Jcosvj; .+k co ( 1 -cos 6 .cosp Jcosif; . 
Jr*xj s s J Jrxj j •t-’yj ss J J jss jj j 

( 14 . 4 ) 

f® .-k w^sinib .=f^ .cosd; .+{ff .-k w^cos 6 .cosS ,)sindj ,-k w^(l-cos 6 .cos 6 .)sirul/ . 
Py 3 s s ^3 J Pyj s s 3 '^ 3 ' ^3 s s' 3 ^3 J 


f .= £Z . 
Pzj Pzj 


(14. 5) 
(14.6) 
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I 2 

Note in Eqs. (14. 4) and (14. 5) that the large k w terms have 

s s 

been subtracted from the left side of each equation and that the remaining 

I 2 

k CO terms on the right side are each multiplied by (1-cos 6 .cos^ .), 
s s 3 3 

which is small compared with unity. 

Next we sum the dimensionless force components due to the jth 

blade, as given in Eqs. (14.4), (14.5), and (14.6), over the N blades. 

We continue to use 2p m ^ as the normalizing factor. Thus 

c s mcix p max 

let 



2p 



m Q ^ 
c s max p max 



2p 



m 

c s max 


pmax 


Rz 2p 



m 

c s max 


pmax 


Then we have 


(14.7) 



, N ^ k CO N 

= — S (f^.+k CO cosijj .)- ^ S cosi]; . 

N j=l ' Pxj s s N 


(14. 8) 


f® - N k' N 

= — S (f® .-k* CO sini|j .)+ Zy simj) . 

N N ''pyj s s N 


(14.9) 


f® N 

Rz 1 •y' 

N “ N . , 

J=1 




ZJ 


(14. 10) 


But because the N blades are equally spaced in azimuth along 2it radians, 
■N N 

s cos»1j .=0 and Zy sin4< .=0. I. e. , tixe principle centrifugal force terms 

j=l J j=l J 
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acting on each blade cancel out* Only if a blade is lost do these terms 
need to be included in Eqs. (14* 8) and (14* 9). Thus optimal scaling is 


preserved* 

Finally, the S-frame dimensionless force components in Eqs. 

(14. 8), (14* 9), and (14. 10) above are transformed to hub and hence body- 

T 

axis components using the matrix [ T^^] , where [ T^^] is given in 

Eq, (2, 3). Thus 


1 i 1 -S . , 1 r® ■ - - 

- Tt f-n cos 1 - + — f_ sm 1 . sm x. +7;f„ cos sin i. 

N Rx N Rx 9 N Ry 0 NRz 0 


(14. 11) 


13 1 .s . 1 _s . . 

Tt f cos I - r; f ^ sin i . 

N Ry N Ry cj> N Rz cj> 


(14. 12) 


1 -b 1 -S . . , 1 .S . . ■ , 1 r® - • /ivi 

“ f = - r; f _ sm i^+r; f^^ sini ,cos + — f „ cos i, cos i. (14* 13) 

N Rz N Rx 0 N Ry <|) 0 N Rz (J) 0 ' • / 

B.15. Block 15 

The hinge force components . and F^ . along the x . and z , 

^ Pxj pzj ^ rj rj 

R-frame axes for the jth blade were normalized in Section B.14 by dividing 

by 2p m ^2 ^2 , This suggests that we normalize the jth blade 

^ c s max p max 

T r 

hinge moment components M, . and M, . by dividing by 2ep m ^2 ^2 

^ nxj nz j ^ ® ^ c s max p max 

Thus let 


M, 


^ hxj 

nxj 2ep ml2 S2 


M. 


^zj = T; 


hzj 


p m ^2 ^2 

c s max p max 


(15,1) 


c s max p max 

Then from Eqs* (15. 2) and (15. 3) in Section A. 2, 15 the equations for 


dimensionless moment components at the jth blade hinge become 
r 

”^hxj ■ f M ^ i ■ M 


■A K, n p. 

^ P. + (15.2) 


pmax 


f2 

p max s max 
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(15. 3) 


r 






M 


pmax 


[ f(5 j)+ sgn 6jF g ] - ^ 


xj 


where 


M 

pmax 

h ~~ 2ep m ^ ^2 

^ c s max p max 


(15-4) 


Here the terms inside the square brackets on the right side of Eqs, (15.2) 

and (15. 3) have already been computed in Eqs. (11. 1) and (11. 2), Section 13. 11. 

The dimensionless moment components in Eqs. (15. 2) and (15, 3) 

T 

are resolved into S-frame axes using [ » where [ is given 

Y 

in Eqs. (2.10). Thus, since we have 


TTL . = m- . sinib . 
hxj hxj ^ j 

s r , 

m, . = m. .cosip . 
ayj 3 

s r 

m_ . = m, . 
hzj hzj 

Next we sum the hub-moment components of Eqs. (15. 5), (15.6), and 
(15. 7) over all N blades to obtain the total dimensionless hub moment 
components. Thus 

s 


(15. 5) 

(15.6) 

(15.7) 


N 


_ JL ^ s 
- N ""hxj 


(15. 8) 


m 


^ = 


N 


N 

1 V s 


s 

_ 1 

N 

E 

N 

N 

j=i 

m ^2 




s 

“h 


ZJ 


(15.9) 


(15. 10) 


s max p max 


is still used as the normalizer 
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for the dimensionless moment components m , m , m . These are 

nx ny nz 

resolved into hub and hence body axis components using the matrix 


[ > where [ is given in Eqs. (2, 3). Thus 


sh 


b s s s 

m, m^^ 

ny hz 

cos i^ + sin i^ sin i^ + — cos i^sin i^ (15.11) 


N 


N 


0 N 


b s s 

m, m. Tx\ 

hy hy . _ hz 

‘ cos 1 , - ‘“TT” sm 1 , 

N N N 


(15.12) 


o 

N 


s 

”hx . . , 

~ ‘e 


m. 


hi 


m. 


N 


sin i, cos i^ + 
<)) 0 


hz 


N 


cos i , cos i^ (15. 13) 
cp 0 


The total rotor moment about the aircraft c, g, is the sum of the hub 
moment plus the moment due to the rotor force acting at the hub, i. e. , as 
given in dimensional form in Eqs. (15.16), (15.17), (15.18), Section A. 2. 15. 

To normalize the total rotor moment components about the air- 
craft c. g. we use 2Nd p m ^ ^2 , where in Section B. 4 we 

max ^ c s max p max 

defined s.s the largest of the relative hub to c. g. position -vector 

components, Hence we let 


b 

Rx 


Rx 

2Nd p m^2 ^2 

max c s max p max 




2Nd p m^2 ^2 

max^ c s max pmax 


b 




z 



2Nd p ml2 Q 

max^ c s max p max 


(15. 14) 


Then Eqs. (15.16), (15.17), and (15.18) in Section A. 2. 15 become 
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= (“V 
Rx ' d 


b / 

“^Ry = 


,b b 

Rx ^Rz 

N "‘d ^ N 

max 


b e , “ha I / , ’'h 

“Ra- >“+'— > N -<T“> 

max pciax max 


(15.15) 


(15. 16) 


(15.17) 




\z 

\i 

^hx’\y’\z 

A . 
ri 

aS aS aS 

A , A , A 
rx ry rz 

a’^ .,A^ .,A^ . 
rxj' ryj' rzj 

A . 

..a’^ .,a^ . 

pxj’ pyj pzj 

aP ..aP .,aP . 

pxj pyj’ pzj 


Cl. 

aP ^aP .,aP . 

cxj cyj czj 


^IS 

%x* %y'%z 

a - a .a 
hx* hy* hz 

s s s 
a , a , a 
rx ry rz 

r r r 
a .,a .,a 

rxj' ryj rzj 

r r r 
a .&a .^a . 

pxj pyj pzj 
aP .,aP .,aP . 

pxj pyj pzj 

aP .,aP .,aP . 
CXJ cyj’ CZJ 

B -frame 


IS 


LIST OF SYMBOLS 

Acceleration of B-frame w, r.t. I-frame 
Components of A along x, ,y, ,z, body axes. 

Dl D D D 

Acceleration of H-frame w. r. t, I-frame. 

Components of A^, along 

Acceleration of R -frame w. r.t. I-frame. 


Components of A . along x ,y , z axes. 

ri s s s 

Components of A . along x - y z . axes. 

ri ® rj, rj rj 

Acceleration of P. -frame w. r.t, I-frame. 

J 

Components of A . along x . y .,z . axes. 

Pij rj,^rj rj 

Components of A . along x ,,y z . axes. 

^ pij ^ PJ P3 

Acceleration of the c. g. of the jth blade w, r. t. I-frame. 

Components of A . along x ,,y z . axes. 

cij PJ PJ PJ 

Constant in the blade geometric pitch equation. 

Dimensionless components of A along x, y, z axes. 

hj 

Dimensionless components of A^. along 

Dimensionless components of A . along x ,y ,z axes. 

rl s s s 

Dimensionless components of A . along x . y .,z .axes. 

ri rj, rj rj 

Dimensionless components of A . along x .,y .,z . cixes. 

pi. ® rj 

Dimensionless components of A . along x .,y .,z . axes. 

Pij PJ PJ PJ 

Dimensionless components of A . along x .,y z . axes, 

cij ^ PJ ^PJ PJ 

The aircraft body-axis frame (x^, y.^, z^) . 

Constant in the blade geometric pitch equation. 
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^djk 

Sjk 

°jk 

d 

max 


F F F 
bx* by' bz 


F^ 

Rx ’ Ry’ Rz 

IT® IT® IT® 

^Rx’ ^Ry’ ^Rz 
^Tjk' ^Rjk’ ^Pjk 


Pj 

^Pxj' ^Pyj’ ^Pzj 

F’^ F’^ 

■*' Pxj’ Pyj' ■*' Pzj 

F® F® F® 
^Pxj’ ^Pyj' ^Pzj 

^x' ^y’ ^z 
b fb 

Rx’ Ry’ Rz 

fS fS 

^Rx’ Ry’ Rz 
^Tjk’ ^Rjk’ ^Pjk 


^Pxj’^Pyj’^Pzj 

if f ^ f ^ 
Pzj 


Drag coefficient for kth segment of jth blade. 

Lift coefficient for kth segment of jth blade. 

Drag of the kth segment of jth blade. 

Largest of x , y , z (hub displacement from c. g. ). 
irl Jri rl 

Position of hinge point for jth blade w. r. t. the hub. 
Magnitude of e^. 

Total external force acting on aircraft, less gravity. 
Components of F^ along ^ body axes. 

Total force on aircraft due to rotor system. 

Components of F^ along y|3» ^ body axes. 

Components of F along x ,y ,z axes. 

XV S S S 

Aerodynamic forces for j blade, kth segment along 

X . , y . , z . axe s . 

PJ PJ PJ 

Hinge shear force. 

Components of F_ along x .,y .,z . axes. 

Pj PJ PJ PJ 

Components of F along x .,y .,z . axes. 

rj rj rj 

Components of F along x ,y ,z axes. 

, S B S 

3 _ 

Dimensionless components of F along x , y , z axes. 

b b b b 

Dimensionless components of F along x, , y, , z, axes. 

XV b b b 

Dimensionless components of F^^ along x^,y^z^ axes. 

Dimensionless components of aerodynamic force for 
jth blade, kth segment, along x^, y^, z^ axes. 

Dimensionless components of F_ along x .,y .,z . axes. 

Pj PJ 'pj pj 

Dimensionless components of F_ along x .,y .,z . axes. 

^ P^ ° rj ''rj rj 
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f® f® 
^Pxj'^Pyj’^Pzj 


f( 5 .) 


G ,G ,G 
X y z 

G®,G®,G® 

X y z 

G^.,G^.,G'. 
xj yj zj 

gP..gP.,gP. 

xj' yj’ zj 
g„ 


®x‘V®z 

4'gy'®z 

4j*4j’4j 




Dimensionless components of F along x ,y ,z axes. 

Jt^a S S S 

J 

Constant in eqiiation for hinge lagging moment. 

Function in equation for hinge lagging moment. 

Gravity acceleration vector. 

Components of G along Xj^,y^,z^ axes. 

Components of G along x ,y ,z axes. 

s s s 

Components of G along axes. 

Components of G along x .,y ,,z . axes. 

PJ PJ PJ 

Sea level gravity acceleration magnitude. 

Dimensionless components of G along x^, y^, z^^ axes. 

Dimensionless components of G along x ,y ,z axes. 

X s s 

Dimensionless components of Galong x ,,y .,z . axes. 

rj rj rj 

Dimensionless components of Galong x .,y z . axes. 

PJ PJ PJ 

Angular momentum of aircraft about c. g. (not in- 
cluding the rotor blades). 


H . 

P^ 


Angular momentum of jth blade about hinge point P 


J 


H -frame 
h 

I-frame 


The hub -axis frame, 

Altitude 

The inertial reference frame (non accelerating 
atmosphere). 


[ij 


Aircraft inertia matrix about c. g. (not including the 
rotor blades). 


I ,I ,I 

XX yy zz 






Principle moment of inertia components of [ I] along 

Blade inertia matrix about hinge. 

Components of [ I^] along x^, z^ axes. 


102 


I 


i. »i .i ..i • 

D s rj pj 




K ,K 
“l "2 

k_ , k , k k . 
b* s* rj' PJ 


Unit vectors along x. ,x ,x .,x . axes. 

s* rj PJ 

Pitch and roll angles of shaft axes w. r.t. hub (and 
hence body) axes. 

Index denoting blade number. 

Unit vectors along y, ,y ,y .,y . axes. 

^ *^b ^rj '^pj 

Constants in the blade geometric pitch equation. 

Constants in equation for hinge flapping moment. 

Unit vectors along z, , z , z z . axes. 

b s rj PJ 

Scaling constant, Eq. (13. 9), Ap. B. 

Scaling constant, Eq. (5. 6), Ap. B. 

Scaling constant, Eq. (13,4), Ap. B. 

Scaling constant, Eq. (1.10), Ap. B. 

Scaling constant, Eq. (9.7), Ap. B. 

Scaling constant, Eq. (1.10), Ap. B. 

Scaling constant, Eq. (15.4), Ap. B. 

Scaling constant, Eq. (1.15), Ap. B. 

Scaling constant, Eq. (9.7), Ap, B. 

Scaling constant, Eq. (7. 7), Ap. B. 

Scaling constant, Eq. (7.7), Ap. B. 

Scaling constant, Eq, (4. 9)t Ap. B. 

Scaling constant, Eq. (8.5), Ap. B. 

Scaling constant, Eq. (13.9), Ap. B. 

Scaling constant, Eq. (11, 3), Ap. B. 

Aircraft mass, not including rotor blades. 
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'•SB'* 



Total external moment acting about aircraft c. g. 



Components of along ^ axes. 



Moment about aircraft c.g. due to rotor system. 



Components of along y^, axes. 



Moment about hub due to rotor system. 


^hx' ^hy' ^hz 

Components of M^ along x^,y^, axes. 


^'^y'"^hz 

Components of M. along x ,y ,z axes. 

n s s s 


^h. 

Moment about hub due to jth blade. 


Klj'Kyj'Kzj 

K^yKyyK^i 

Components of M along x ,y ,z axes. 

Xl. s s s 

_ J 

Components of M. along x ,,y .,z . axes. 
^ h. ^ rj’ -^rj rj 

Flapping hinge moment for jth blade. 


FDj 


T^• 

LDj 

Lagging hinge moment for jth blade. 


M„ 

Jk 

Mach number of jth blade, kth segment. 


M 

aero 

Blade moment about hinge due to aerodynamic forces. 

M, . 
hinge 

Blade moment about hinge due to hinge moments. 


“Bj 

Total moment acting about hinge on jth blade. 


M . 

PJ 

Total moment (excluding gravity) acting about hinge 
on jth blade. 

M . M . 
pxj, pzj 

Components of M_ along x z . axes. 

p PJ PJ 


m 

Blade mass 


"\x*^y*”^bz 

Dimensionless components of M^ along x^, y^, z^ 

axes. 

■"Lc-”Ry-”R 

' z 

Dimensionless components of M^ along x^ , y^ , z^ 

axes. 


Dimensionless components of M^ along x^,y^, z^ 

axes. 

"'L-”hy'“hz 

Dimensionless components of M along x ,y , z 

xl S S S 

axes. 
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’^yj' 

”hrxj' '^yj' ”Hizj 


m m 
pxj - pzj 

N 


P. -frame 
3 

^b'°b’^b 

PJ.qJ.r8 

p^..q'.,r5 . 

rj rj rj 

pP.,qP.,rP. 

rj' rj rj 

P .,Q .,R . 
PJ PJ PJ 

W'b 

s s s 

Pr'V^ 

P^..q^.,r^. 

rj*^rj’ rj 

P P P 
rj 

P .f<l -»r . 
PJ PJ PJ 

R . -frame 
J 


S-frame 


S ,S ,S 
X y z 


t VI 


fVi 


tvl 


Dimensionless components of along x^.y^.z^ axes. 

_ j 

Dimensionless components of along x^^,y^j,z^j axes. 

Dimensionless components of M . along x z . axes. 

PJ PJ PJ 

Total number of blades in rotor system. 

The jth blade span axis frame, * s»yn 4 »^u 4 * 

PJ PJ PJ 

Components of along axes. 

Components of along ^q*Vq 9 axes. 

Components of ^2 . along x .,y z . axes. 

^ ri rj rj rj 

Components of ^2^^ along 

Components of ^2^^ along ^pj»^pj*^pj ^®®- 
3 

Dimensionless components of along x^,y^,z^ axes. 

Dimensionless components of ^2^^ along x^,y^,z^ axes. 

Dimensionless components of along Yj.j» 

Dimensionless components of ^2^. along ^pj'^pj’^pj 

Dimensionless components of along ^pj^^pj^^pj 

J 

The jth rotating shaft-axis frame, x^^, y^j,z^j. 

Position of hub w. r. t. aircraft c. g. 

The shaft axis frame, x , y , z . 

o S S 

Distance north, east, and downward. 

Number of segments per blade. 

Matrix for transformation from body to earth axis 
components. 

Matrix for transformation from earth to body axis 
components. 

Matrix for transformation from H to S-frame axis 
components . 
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f 


Matrix for transformation from S to H-frame axis 
components. 


[ T ] 

^ sr . 

J 


Matrix for transformation from S to R. -frame axis 
components . ^ 


[ T ] 

rs. 

J 


Matrix for transformation from R. to S-frame axis 
components . ^ 


[ T ] 


Matrix for transformation from R. to P. -frame axis 
components. ^ ^ 


[ T ] 

pr. 

J 

U V w 
b’ ^b’ b 

hb’ hb' hb 

U®, V®, W® 
r r r 

r r r 

rj' rj’ rj 

PJ PJ PJ 

U^., V^.,wP. 

PJ PJ PJ 

^Tjk’ ^Pjk' ^Rjk 

V\*^b 

w ^ 

hb* hb* hb 


u 


s „s 


, w 


. s 


rs' rs' rs 

r r r 

u .,v .,w . 

PJ PJ PJ 

P P P 

PJ PJ PJ 

^Tjk'^Rjk’ ^Pjk 


bi 


hb 

hi 


Matrix for transformation from P. to R. -frame axis 
components. ^ ^ 

Components of along x^,y^, axes. 

Components of along x^,y^,z^ axes. 

Components of V . along x ,y , z axes. 

ri s s s 

Components of V . along x .,y .,z . axes. 

ri ^ rj*'^rj* rj 

Components of V . along x .,y .,z . axes. 

pi. ® rj ''rj rj 

J 

Components of V . along x .,y .,z . axes. 

Pij PJ PJ PJ 

Components of V., along x .,y ,,z . axes. 

jk ® PJ ■'pj PJ 

Dimensionless components of V^. along x^,y^,Zj^ axes. 

Dimensionless components of along x^^ , y^ , z^^ ctxes. 

Dimensionless components of V . along x ,y ,z axes. 

rl s s s 

Dimensionless components of V . along x .,y z . axes 
^ pi. ^ rj*'^rj* rj 

J 

Dimensionless components of V . along x .,y z . cixes 

Pij PJ PJ PJ 

Dimensionless components of V., along x .,y z . axes. 

jk ® PJ’^PJ’ PJ 

Velocity of B-frame, x^,y.j^,z^ axes. 

Velocity of H-frame w. r. t. B -frame. 

Velocity of H-frame, x^^, y^, z^ axes. 
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vv=% 


X , y , z 


X ,,y „z . 
rj ' t ) rj 

I t t 

X .,y .,z , 
r J r j rj 


X y .,z . 
P3 PJ PJ 


TRANS., 

Jk 


Velocity of R -frame, x .,y z .axes, 

^ rj ^ r j r j 

Velocity of P. -frame, x .,y .,z . axes. 

3 PJ PJ PJ 

Velocity of c.p. of kth segment, jth blade. 

Magnitude of V^j^ 

Dimensionless V,, 

Jk 

Body axes (B -frame) 

Hub axes (H -frame) 

Shaft axes (S -frame) 

Rotating shaft axes, jth blade (R^ frame). 

Rotating shaft axes with origin translated to 
Blade-span axes for jth blade (P^ frame). 

Displacement from jth blade hinge to kth segment c, p. 
Angle of attack, jth blade, kth segment. 

Angle of attack for drag and lift coef. computation. 
Flapping angle, jth blade. 

Angle of attack yaw angle, jth blade, kth segment. 

Swash plate rotation angle. 

Lagging angle, jth blade. 

Body axis pitch angle. 

Blade twist angle, jth blade, kth segment. 

Yawedblade twist angle, jth blade, kth segment. 
Displacement of jth blade, kth segment c. p. w. r. t. hinge. 
Displacement of jth blade c. g. w. r. t. hinge. 
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■(>. 








Q 


H 



^ . 
SI 


^ . 
ri 



Q 

s 

pmax 


CO 

s 


Body -axis bank angle* 

Aerodynamic angle of attack, jth blade, kth segment. 
Body -axis heading angle. 

Azimuthal angle for jth blade -span axes. 

Angular velocity of B -frame. 

Angular velocity of H-frame (same as B-frame). 

Angular velocity of S-frame (same as B-frame). 

Angular velocity of R. -frame. 

3 

Angular velocity of -frame. 

Angular velocity of R^-frame w. r. t. S-frame. 

Normalizing factor for dimensionless blade -span 
axis angular velocity components. 

Dimensionless ^ . 

s 
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